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I. INTRODUCTION 



A. CFAR SIGNAL PROCESSING 

Signals processed by radars are inherently corrupted by environmental clutter and 
thermal noise from receiver components. The echo signal returned from any potential target 
will thus be accompanied by noise from these sources. The simplest method of detecting 
the target return signal is to set a fixed power amplitude threshold level which is known to 
be greater than the mean noise level and most of the noise peaks, but less than the expected 
magnitude of the target return signal. Thus, any signal which is measured to be greater than 
the threshold is declared by the radar to be a target detection, as shown in Figure 1. Such a 
detection may include a noise peak, and thus a false detection. This is known as a false 
alarm. By measuring the average number of false alarms occurring in a specified period of 
time, a false alarm rate can be established. The ratio of the average number of false alarms 
to the total number of noise samples is known as the probability of false alarm (Pfa)- By 
the same token, the ratio of the average number of target signal detections to the total 
number of target signals sampled is known as the probability of detection (Pj)). This 
implies that some target signals may not be detected, and thus the probability of a missed 
target detection can be computed. In general, the detection performance of a radar is stated 
in terms of the PpA and the P^. 

Fixed threshold detection is inadequate for radar systems that have to detect the 
presence of targets within background environments that are more complex and less known 
than thermal noise alone. This type of detection produces large changes in false alarm rates, 
a very undesirable characteristic for radar systems intended for high target detection 
sensitivity and reliability. Also, the great number of detections that occur saturate and 
disrupt the tracking computer associated with the radar system. Thus detector design 
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criteria for modem military radars is to maintain a low and constant Pp^ as weU as a good 
Pp) for a given signal-to- noise power ratio (SNR). Automatic detection systems which 




threshold 



mean 

noise 

level 



time 




Figure 1: Detection Threshold 



maintain a constant Ppy^ are known as constant false alarm rate (CFAR) systems. In these 

systems, the constant false alarm rate is maintained by adaptively varying the detection 
threshold level in accordance with estimates of the noise power (targets not present). Such 
detectors can easily cope with fluctuating noise levels, as the estimates are taken in real 
time. CFAR detection is usually performed in range. That is, the noise estimates are taken 
from a predetermined number of range bins (also known as cells) which are adjacent to the 
range bin on which detection is being performed (known as the test cell). This is shown in 
Figure 2. These adjacent cells are called reference cells. The analog-to-digital converter 
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samples taken over the time period occurring between two consecutive transmitted pulses 
are the range bins. That is, each range bin represents an increment in range from the radar. 



transmitted radar pulses 




Figure 2: Unambiguous Radar Range with a Sample CFAR Window 



The maximum possible range for the time period between the two pulses is known as the 
unambiguous range. The reference cells and test cell form a window which slides through 
a specified set of range bins, evaluating each range bin for a target. 

Most adaptive threshold detectors assume that the samples in the reference range cells 
are independent and identically distributed. Furthermore, it is usually assumed that the time 
samples are independent. In essence, the detectors test whether the test cell power is 
sufficiently larger than the noise power in the reference cells. It is assumed that the 
probability density function of the noise is known except for a few unknown parameters. 
The suiTounding reference cells are then used to estimate the unknown parameters, and a 
threshold based on the estimated parameters is obtained. Since the determined parameters 
are estimates, the resulting threshold has some error and must therefore be slightly larger 
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than the threshold that one would use if the parameters were known exactly a priori. This 
causes a loss in target detection sensitivity, called a CFAR loss, and thus results in a 
degradation of the Pj). 

In order to analyze the Pd and the PpA performances theoretically, the probability 
density functions of the radar range bins have to be derived. In general, given a density 
function for a non-target cell (noise only) and one for a target cell (signal with noise), the 
PpA and the Pp are obtained by integrating over these functions for a given detection 
threshold T, as shown in Figure 3. Obtaining these functions can be very tedious. The 




Figure 3: Pd and Pp^ 



resulting analytical expressions representing these density functions are not always 
practical to manipulate and program. Thus, attempts are often made to model the equations 
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with simple closed form expressions to facUitate the analysis and implementation of the 
desired CFAR processing scheme. 

B. PRINCIPLE CONTRIBUTIONS 

The main contribution of this research is the rigorous treatment of the effects of two 
envelope detection approximations on the statistical characteristics of the detected signals 
within a greatest-of CFAR (GO CFAR) signal processor. The first approximation is given 
by Xg = aMAX { |/|,| Gl } + bMIN { |/U (21 } . where I and Q are, respectively, the in-phase 
and quadrature samples of the input signal. The values a and b are simple multiplying 
constants. A closed form solution for the probability density function of this detector is 
shown to be untrackable analytically. An alternative method to obtain a closed form 
expression for the density function is developed. This approach curve-fits the results of a 
numerical convolution. These closed form expressions enable a quick implementation of 
the function within system simulation programs. Reference tables provide the necessary 
data for this method. 

The second approximation investigated is given by Xg - a|/l +h|(2l, a simplified 
form of An analytical closed form solution is presented, as well as a curve-fit solution 
similar to that obtained for Xg. This allows a comparison of the two envelope detection 
approximations in terms of their coefficients. 

Finally, a comparison between the true envelope detector GO CTAR performance and 
the performance using x^ and kg is made in terms of numerically derived probability 
density functions and the Pj). Results show that there is a significant difference in 
processor performance for the case of I and Q being zero mean white Gaussian noise with 
a non-fluctuating target. 

C. THESIS OUTLINE 

The objective of the research documented here is to investigate a method of 
computing the Pj) analytically for a GO CFAR processor that uses the envelope detection 
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approximations = aMax {\I\,\Q\} + bMin {\I\,\Q\) and = a|/| + Z?I2I • The main 
hurdle in doing this is to determine a closed form expression for the probability density 
function of a detected radar range cell with a target signal present in it. 

Section II describes the operation of a GO CFAR in a simple manner. The 
functionality of it is described by stepping through all its operations sequentially, providing 
a brief discussion at each point. It then provides a short discussion of envelope detection 
approximation. Finally, the requirements for the analytical determination of the Pd for the 
GO CFAR detector is outlined. Section m provides the rigorous analytical development of 
the closed form solution for the probability density function of a target test cell using the 
envelope detection approximation x^. It is a step-by-step walk-through of the derivation, 
presented first as a general result and then applied to the particular requirements of this 
research. The difficulties encountered and the assumptions made for every step of the 
derivation are described. This section also provides a closed form expression derived from 
curve-fitting the numerical convolution. Section IV presents the closed form solution for 
the probability density function of a target test cell using the envelope detection 
approximation Xg. A closed form solution using curve-fitting techniques is also described. 
Section V presents a comparison of the two methods of envelope detection approximation 
discussed in Sections HI and IV in terms of the Pj) and the PpA- Section VI provides 
concluding remarks on the success of this research. It also discusses the possibilities for 
continued research in this field, as well as different aspects of the problem which could be 
further studied. 

D. REFERENCE OVERVIEW 

The studies by P.E. Pace and L.L. Taylor [Ref. 1] concentrate on obtaining the PpA 
analytically for an envelope detection approximation GO CFAR processor. The 
approximation used is Xg = a\I\+b\Q\. They derive a closed form solution of the 
probability density function of a radar range cell with noise only in it. Being able to 
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integrate this function in closed form, they are able to produce Pp^ curves numerically. The 
curves are derived as a function of a threshold multiplier and the number of reference 
cells used on each side of the test cell. They also derive, through curve-fitting, closed form 
expressions for the Pp^ and compare its performance to the numerical solution. 

The technical memorandum by D.J. Wilson [Ref. 2] derives the Pp for a test cell with 
an envelope detection approximation with a = 1 and h = 1. He provides a useful insight 
into the analysis of this type of approximation and the subsequent derivation of the Pp. 

Another technical memorandum, this time by D.G. Loberger [Ref. 3], compares the 
performance of a cell-averaging CFAR processor (similar in concept to GO CFAR 
processing) for two different configurations. He analyzes the situation in which the test 
range ceU is not used in the noise averaging of the reference cells. He also examines the 
situation in which the test cell is included in the averaging. The important result of this 
paper is to relate the threshold multipliers of each situation considered with a simple 
conversion formula. Using this result as a constraint guarantees the same Pp^ and Pp, 
performances for either of the envelope detection approximations. 

A study by G.P. Laulusa [Ref. 4] evaluates the PpA and Pj) performances of a GO 
CFAR processor which includes the test cell along with the reference cells in the estimation 
of the noise level. He uses the envelope detection approximation with unity multiplying 
coefficients and a normalized version of the probability density fimction of the test cell. 

The paper by V.G. Hansen and J.H. Sawyers [Ref. 5] determines the CFAR loss of a 
GO CFAR scheme as compared to a simple cell-averaging CFAR. In the process of doing 
this they derive an integral expression for the Pp in terms of the probability density function 
of a greatest-of output block. This, in effect, defines the main building block for any Pp and 
PpA computations on a GO CFAR processor, for any type of signal detector at the input of 
the processor and for any variation of the GO CFAR processing scheme. 
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II. GO CFAR OVERVIEW 



A. IMPLEMENTATION 

A typical GO CFAR processor is shown in Figure 4. It uses an envelope detector, a 
lead window of reference cells, a lag window of reference cells, a test ceU, two adders, a 
maximum value detector, two multipliers and a binary comparator. The true envelope 

detector is x = Jp' + q} , where / and Q represent the in-phase and quadrature samples, 
respectively, of the input signal. . The I and Q samples are obtained from an analog-to- 
digital converter (ADC), which produces the digital representations of the radar return 
signal. The ADC receives its analog I and Q signals from the synchronous detection and 
down conversion circuits. 

The envelope detected signal samples (from the I and Q samples) are not actually 
written into a GO CFAR processor per say, but rather are stored in consecutive range cells 
(or bins). This is not depicted in Figure 4. The GO CFAR processor is a sliding window 
which reads in the signal samples of the necessary range cells for processing. The sliding 
window enables the GO CFAR processor to test each cell for potential targets. It is divided 
into lead and lag segments, each of length n cells. Selection of n depends on the desired 
quality of the noise estimate, and usually lies between 4 and 32. 

The signal samples of the lead and lag windows are summed and the two results 
compared. The maximum sum is selected to calculate the estimate of the noise level. This 
value is then multiplied by a factor of 1/n to provide the final time-average noise estimate 
for the signal in the test cell. The estimate is then multiplied by T, a threshold multiplier, to 
establish a final detection threshold Vg. This multiplier is preset to a value which will yield 
a desired, constant Pp^. Finally, the detection threshold is compared to the signal present 
in the test cell. If the signal in the test ceil is larger than the detection threshold, a target 



8 



Test 

CeU 



Increasing Range 



► 




Figure 4: GO CFAR Processor 
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detection is declared. It can be concluded that GO CFAR processing is fairly simple in 
principle. However, it is much more complex to analyze theoretically. 

B. ENVELOPE DETECTION APPROXIMATION 

The envelope detector can be implemented with the envelope detection 
approximation mentioned in section I of the form 

x^ = aMax{\I\,\Q\}-\-bMin{\I\,\Q\} . (l) 



Another similar approximation is 



x^ = a\I\ +b\Q \ . 



( 2 ) 



The constants a and b of the envelope detection approximation are predetermined so as to 
yield a desired performance level for the approximation given in (1). Envelope detection 
approximation is performed mainly for two reasons. First, much less hardware is required 
for implemention (less expensive). The squaring operation and the square-root operation 
in the true envelope detector are also difficult to design digitally so as to yield answers with 
minimal numerical rounding errors. 

Secondly, envelope detection approximations perform faster in real time compared to 
the true envelope detector. This is a crucial requirement, since the sequential GO CFAR 
processing on all the range cells must be completed prior to the next received pulse. For this 
reason, GO CFAR processors are frequently found in low pulse repetition frequency 
(LPRF) radar modes of operation. 

The Max and Min operators in the envelope detection approximation x^ have the 
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effect of uniformly distributing each complex input sample, of the form 5 = I +jQ, over 
the range of angles 0 = 0... (ti/ 4) radians. This dependence on 0 can be expressed by 
rewriting the Max and Min terms in polar coordinates as 

X = x^(a • cosd + b • sind) (3) 



The relative approximation error can be defmed as 

e(0) = — = 1 - (a • COS0 + Z) • sin0) (4) 

Then, the performance criteria in terms of the average error and mean-square error can be 
expressed as, respectively. 



n/A 

e = i|e(e)d6 (5) 

0 



^ J e^(Q)dQ . (6) 

0 

A number of envelope detection approximations are shown in Table 1. [Ref. 1]. 

The coefficients of Table 1 are considered in this research. This allows the findings 
fi-om this thesis to be correlated to the results of previous research in this area [Ref. 1]. They 
provide a good variety of performances, in that the errors can be observed to vary over a 
wide range. The sixth approximation is such that it yields a zero average error. 
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TABLE 1: ENVELOPE DETECTION APPROXIMATIONS 



approxi- 

mation 


a 


b 


e 




1 


1 


1 


-213 


30.0 


2 


1 


1/2 


-8.68 


9.21 


3 


1 


1/4 


0.65 


4.15 


4 


1 


3/8 


-4.02 


4.76 


5 


31/32 


3/8 


-1.2 


2.7 


6 


0.948 


0.393 


0 


2.33 


7 


0.96043 


0.39782 


-1.3 


2.7 
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C. ANALYTICAL DETERMINATION OF PROBABILITY OF 
DETECTION 

It has been assumed for this paper that the 7 and Q input samples to the GO CFAR 
processor are independent and identically distributed. It is also assumed that the noise 
components of the samples are normally distributed zero mean and unity variance 
(~N(0,1)). Once the probability density function of the target test cell has been obtained, 
the probability of detection Pp can be determined using 



Pd{SNR) 




dz 



(7) 



The probability of false alarm Ppy^ is given by [Ref.l] 
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n 



( 8 ) 



The term ^ represents the density function for the n reference cell summation and 
represents the density function of the test cell. From Figure 4, ^ is obtained by 

performing an n-fold convolution of the density function of a single reference cell. That is, 
each reference cell is assumed to have the same noise distribution and to be independent of 
each other. The distribution of the term p^ is the only difference between the Pp (SNR ^ 0, 
target present) and the PpA (SNR = 0, no target, noise only). In the former case it is 
dependent on SNR, whereas in the latter case it is not. 

The true envelope detector would yield a Rician distribution (or non-central 
Rayleigh) for the target signal in each of the range cells. Closed from expressions for the 
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probability density function (pdf) of such a distribution exist and involve modified Bessel 
functions [Ref. 8]. The true envelope detector yields a Rayleigh distribution for noise alone. 
However, since the envelope detection approximations and Xg are being used, the 
derivation of the pdf turns out to be much more difficult. This forms the focus of this 
research. 

The goal is to attempt to determine if a closed form solution can be derived for the 
pdf of the signal produced by (1) and (2). The approach to be taken consists of deriving the 
probability density functions of the smallest components of the equation first. Then the 
transformed density functions have to be derived following the same constructions that 
make up (1) and (2). The pdf of Xg is Pg {x) , representing the pdf of a range cell with a 

target signal present in it. The pdf of Xg is Pg (Jc) . It is assumed that the target signal is 
located within the test cell only, and does not contribute power to the noise-only cells 
adjacent to it. 
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III. PDF FOR THE ENVELOPE DETECTION 
APPROXIMATION aMax{l/l,l!2l} + bMin{\I\,\Q\} OF A 

TARGET TEST CELL 



The development outlined in this section derives the density function {x) of a 
target test cell sequentially. It begins with the probability density functions of each of the 
I and Q sample expressions. These are known to be Gaussian distributed, being composed 
of known means and normal noise. The densities of the magnitudes of I and Q are found 
next. Then the density functions of the terms Max {\I\,\Q\} and Min {\I\,\Q\} are 
derived. This is followed by the derivation of the density functions of these terms 
multiplied by constants, i.e. aMax {\I\,\Q\} and bM in {\1\,\Q\} • Finally, the probability 
density function of the term x^ = aMax {\I\,\Q\} +bMin {\I\,\Q\} is derived. The 
derivation of (x) involves a convolution of the two terms in x^. Because the density 

functions of Max and Min operations are quite involved, the convolution of their terms 
generate a very complex probability density function. This may not be the most desirable 
approach. As shown, this density function is such that computational efficiency and 
computer memory capacity become important concerns. Additionally, in order to obtain a 
closed form solution of the convolution integral, curve-fitting techniques are used to model 
the error function. Instead of using what would turn out to be an extremely complex closed 
foim solution for p^ (x) , curve-fitting techniques are also investigated to generate 
expressions for p^ (x) . The expressions for p^ (x) are derived for the selections of a and 
b given in Table 1 . The coefficients necessary to implement these expressions are tabulated 
to provide a useful design tool for the evaluation of the Pp and to provide a performance 

comparison with other envelope approximations such as x^. 
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A. CLOSED FORM PDF BY ANALYTICAL METHODS 



The inputs into the envelope detection approximation block of Figure 4 are the in- 
phase and quadrature components of the radar return signal. The data sample for each 
channel can be expressed as the sum of a modulated amplitude and some noise, as 

/ = AcosQ + rij (9) 



and 



Q - i4sin0 + nQ 



( 10 ) 



where A is the time-invariant signal amplitude and 0 is the signal phase. The noise 
variables, n/ and hq, are assumed to be mutually independent and identically distributed 
Gaussian random variables with zero mean and unity variance (~N(0,1)). By definition, the 
/ and Q terms are therefore also normally distributed random variables which are 
independent of each other and identically distributed. As discussed in section n, the signal 
phase is uniformly distributed over the range 0 ...(k/ 4). The final result for the 
probability of detection is averaged over the phase. That is, the Pj) is calculated as 



n/4 

F,(SNR)=lj 





(x, SNR, 0) dx 



dz 



dQ 



. ( 11 ) 



The average values of I and Q are, respectively, 

nij = /icos0 (12) 
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and 



niQ = /4sin0 



(13) 



That is, the noise samples have zero mean. The signal amplitude is also related to the signal- 

to-noise ratio (SNR) as A = Jlo^SNR . 

1. PDF of a Gaussian Random Variable 

The probability density functions of / and Q are Gaussian and are expressed as 



P[iO 






1 




2<r 



(14) 



and 



PQiq) = 







(15) 



Since the variance is imity, the terms can be omitted in the following derivations 
without any effect. 

2. PDF of the Absolute Value of a Gaussian Random Variable 

Considering only the / random variable, the cumulative distribution function 
(cdf) of / is the integration of the probability density fimction and can be expressed as 

Pj(i) = Pr{I^i} 
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or 






Pj{i)= Je' 



du 



(16) 



which is the probability that a given realization of the random variable / will not exceed a 
threshold i. The cdf of 1/1 is then 



= Pr{\I\^i} 



or 




(17) 



It is known that the probability density function of I/I is obtained by 
differentiating the cumulative distribution function of I/I. Since the integral operation is 
linear, (17) can be expanded into two parts as 

i 

7*1/1 (0 = jpj(u)du 



or 



0 i 

P\i\(i)= jP!(u)du + jpj(u)du . (18) 

-i 0 
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Differentiating this expansion with respect to i, the density function of 1/1 can be expressed 



as 



= f.jPi(u)du + ^.jpj(u)du 
-i 0 



(19) 



and becomes 



P\ilii) = Pi(i) +Pji-i) 



( 20 ) 



The final expression for the pdf of III is then given by 



P\l\ (0 





u(i) 



or 



P\i\(0 = 




( 21 ) 



where u(i) is the unit step function. This has to be present to ensure that the variable i, 
representing values of 1/1, remains positive in accordance with the absolute value operation. 
In a similar fashion , the probability density function of \Q\ can be expressed as 



P\qM) = 



7^ 






+ e 



'-(q + mQ) 



2 -| ^ 



u{q) 



( 22 ) 
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These functions are valid probability density functions since they satisfy the appropriate 
requirements: they are positive for any value of i and q, and they both integrate to unity. 

3. PDFofMax{l/l,l0l} 

To begin, a random variable Z is defined which takes the maximum of two other 
mutually independent random variables X and Y, 

Z = MAX (X, Y) (23) 



The variable Z has a cdf defined in the following manner 

= Pr{Z^z} 



= Pr{MAX(X,Y) ^z] 



or 



p 2 (z)= Pr{X^znY^z} (24) 

In other words, the cumulative distribution function of Z is the probability that both X ^ z 
and T ^ z for any given value z. Because X and Y are independent from each other, 

F^(z) = Pr{X^z} Pr{Y^z} 



or 



Fz(z)= F^iz)FY{z) 



(25) 
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f 

That is, the joint cdf of the two independent random variables X and Y is equal to the 
product of their respective distribution functions, (z) • Fy (z) . 

In general, the derivative of a function which is composed of the product of two 
other functions (w and v, for example) which in turn each depend on the same variable (x) 
can be written as [Ref. 9] 




u 



dv du 
dx ^ dx 



(26) 



Applying (26) to (25) to obtain the density function of Z, and knowing that differentiating 
a distribution function with respect to its argument yields the corresponding density 
function, (z) is obtained as 

jf^U) =fz(2) =|[f^(z)f,(z)] 









or 



/z(2)= /x(z)^y(^) +/y(z)F;^(z) (27) 



This expression represents the pdf of a maximum of two random variables. It can also be 



21 



derived using [Ref. 8] 



/z(z) = ^ ^fix,y)dxdy (28) 



where the term AD 2 is defined as the region of the jc-j plane where z<Z^z + dz, and the 
function f(x,y) is the joint density function defined as 



f(x,y) =/xW/y(3') (29) 

because of the fact that the two random variables are independent. This approach, however, 
is slightly more complicated. 

In the context of this thesis research, the random variable Z is defmed as 

Z = Max(\I\,\Q\) . (30) 

In order to find its pdf using (27), the cumulative distribution functions of the I/I and \Q\ 
random variables are first derived. By definition, these are obtained by integrating their 
respective density functions over the range to z. Using (21), the cdf of p\j\(i) is given by 

Z 

P\i\(z) = jp^j\(i)di 




di 
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1 









+ e 



•di 



or 




The function erf is the error function defined as 



erf{ 



x — m 



o 



) 




(32) 



Similarly for \Q\, using (22), the cumulative distribution function for is given by 



/’la (-’) = 5 




(33) 



Using (27), the density function of Z is then 



= P\I\ (2)^iqi (2) +P|Gi ( 2 ) ^1/1 (z) 



or 






2j2n 



■-(Z-/7I; 



;)^ r-(Z + /n/)2-j'| 

J + ,L 2 J 



{erfi 



z — m, 



') +erf{- 



z + m 



Q 



)} + 
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u(z) 



f |--(z-mg)^ 


1 p-(z + mp) 


1l 


L 2 


J L 2 


J 


< e 


+ e 











z-m, z + m, 

{erf(-Y^) +erf(-^^) } 



(34) 

That is, the probability density function of the maximum of the two random variables I/I and 
\Q\ is a combination of exponentials and error functions. 

4. PDF of Minim, IQI} 

The density function of a minimum of two mutually independent random 
variables is found in a similar manner as that for a maximum, discussed in the previous 
subsection. The random variable W is assigned to be 

W = Min (X, Y) (35) 



where X and Y are two independent random variables. The cumulative distribution function 
of W is expressed as 



F^(w) = Pr{W^w} 



= Pr{MIN{X,Y) ^w} 



or 






( 36 ) 



24 



This means that the cdf of the minimum of X and Y is equal to the probability that either 
X^w or Y^w for any given value w of W. Thus the distribution function can be expanded 
as per the definition of an or operator as follows 

F^(w) = Pr{X^w} +Pr{Y^w] -Pr{X^wnY^w} . (37) 

Because the random variables X and Y are independent, the cumulative distribution 
function becomes 

F^(w) = Pr{X^w} +/>r{y^w} - Pr w} Pr w} 



or 



F^(w)= F^(w) +Fy(w) -F;f(w)Fj,(w) (38) 

which is expressed in terms of the distribution functions of each of the random variables X 
and Y. Differentiating the distribution function of the random variable W with respect to w, 
and thereby applying equation (26) to the third term, the probability density function of W 
can be written as follows 






= +Fj,(w) -F;^.(w)Fj,(h-)] 
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= fxM +fy(^) - [^x(^‘’)/y(^) +^r(H’)/x(^)] 



= L/x(^) -/x(^)^y(’^)] + l/y(^) -/y(>^)^x(^)] 



or 



/m^(»v)= / x(w) [1-F},(w)] +/y(w) [1-F;f(>v)] (39) 



In this thesis it is desired to find the density function of the random variable 



W = Min(\lUQ\) 



(40) 



In accordance with (39), the density function of W is given by 



Pw^'^) =P|7l(w) [l-^lQl(w’)] +/>K2 i(h') [1--P|;i(vv)] (41) 



Using (31) and (33) for the cumulative distribution functions of the random variables I/I and 
IQI respectively, and (21) and (22) for the density fimctions of I/I and \Q\ respectively, the 
probability density function of W can then be expressed as 
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= 



1 




p-(w>-/n,)2-| 




L 2 J ^ 


L 2 J 




w-mQ\ 

f 2 ) 




w + niQ^ 

J 2 > 



} + 



p-(*v-/ng)^ 




--(w + rtig)^ 


L 2 _ 




L 2 J 



{ 1 - ier/i 



IV - m A 

J 2 . 




-it} 



u (w) 



or 



Pw(^) = 



2jlii 



^L— 2 — J 



fW — mrt\ fw + mrt\ 

\ {2-erf\—^\-erf ^ 



J 2 ) \ J 2 



} + 



/ 



p-(*v-mg)^ p-(w'+/7ie)^ 

.L— 



,2.erf[-^yerf[-^ 



u(w) 



( 42 ) 



This defines the probability density function of the minimum of the random variables I/I and 
IQl, expressed in terms of exponentials and error functions. Note that the unit step fimction 
prevents from being negative. 
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5. PDF of aA/fljc{l/l,l0l} and bMin{\I\,\Q\\ 

Generally, if the random variable F is a linear transformation of one other 
random variable X of the form 

Y = aX (43) 

where a is a constant, the probability density function of Y is then given by [Ref. 8] 

fy^y) = W 

where fx is the density function of X. Let the random variable Z' be given by 

T = aZ = aMax {\I\,\Q\} (45) 



The density fimction ( ^) of Z is given by (34). Assume the constant a is always positive. 

Therefore, using (44), the density function of Z' is then 

or 



Pt (^) = 










2 




2 



• erf 


— — UIq 

a ^ 


+ erf 


a ^ 




2 


2 
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z , ^ 




r ,2 








-(- + /ne) 


^ € 


2 


+ e 


2 



■< erf 


a ^ 


+ erf 


a ^ 


> 


2 


2 






( J 




^ J 





u{z) 



(46) 



In a similar fashion, let W be given by 

W = bW = bMin{\l\AQ\) (47) 

where h is a constant. The density function p^, (w) of W is given by equation (42). Assuming 
that the constant b is always positive, and again using (44), the probability density function 
of is given by 

PwM = iPwC^) 

or 

Py^.{w) = 




r ,w 




~ .w . ^ 








2 




2 



e + e ► 



2-erf 



h-V 



f2 



/W 

‘ -r +tn, 



-erf 



f2 



) + 



/ J 
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r 




r 
















2 




2 




1 ^ 




+ e 







2 - erf 



/-W \ 



J2 



-erf 



/W X 

f - + m/ 



J2 



u (w) 



(48) 



6. PDF of aMax{\I\,\Q\) + bMin{\I\,\Q\} 

The envelope detection approximation is also a random variable and is given by 

X, = aMAX{\Il\Q\) ■VbMIN{\I\,\Q\} 



or 

Z' + W (49) 

The density functions of the random variables Z' and W are known and given by (46) and 
(48), respectively. Since I and Q are assumed to be independent random variables, then Z' 
and W are also independent. This is known because of the fact that at any given time, if Z' 
is a function of l/l, then W is automatically a function of \Q\ only, and conversely if Z' is a 
function of \Q\, then IF* is a function of l/l only, by definition of the Max and Min operations 

mxg. 

If two random variables are independent, then the density of their sum equals the 
convolution of their densities [Ref. 8]. Hence the density function ofXg is the convolution of 

the density functions Z' and W . This represents the density fxmction (x) of a target test 
cell. The pdf of x^ is thus given by 

pAx) = Pz- ( x) (x) 
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or 



PcW= \Pziy)Pw'^^~y)^y 



This linear convolution integral must be solved to obtain the pdf of x^. Since P 2 ' (“) = 0 
for z < 0 and (w) = 0 for w < 0, then Pf.(x) = 0 for jc < 0. In addition to removing 

the need for the unit step function in the density function expressions of Z' and W , the 
integration limits can be set from 0 to a: . The density function of Xg then becomes 



Pc W 




■ erf 


-m ^ 

1- ^Q\ 


+ erf 


a ^ 




2 


2 






^ ) 




^ J 





” V ^ 




r V 'h 








-(^+/ne) 




2 




2 






+ e 


L J 





- erf 


a ' 


+ erf 


(-+">A 

a ' 




2 


2 






V ) 




1 J 





X 



b2j^ 



r i^~y 






2 


+ e 
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2-erf 



rx-y ^ 
h~^Q 



J2 



-erf 



r^~y , ^ 



J2 



+ 



r 








< e 


2 


+ e 


2 



X 



2 -erf 



r^~y ^ 

S--^I 



72 




■dy . 



( 51 ) 



As demonstrated, this probability density function is a complicated expression. 
An evaluation of this expression is presented next. Each term within the integral is 
expanded. This expansion is necessary to enable a closed form solution to be developed. 
The pdf of Xg therefore becomes 






J [ (- 2EQ - 2EH + 2BJ + 2BS - 2EO + 2LO - 2KA A2BW + 
0 



KIR + 2LM - 2KW + 2LN - 2KU - 2EG - LOG - BXW -BDW-¥ 2LG - 



LPN + KPS + 2LQ - 2KV - BDA - 2EM + EIO + KIV - LIH + KPU- 
BPS + EXM - LIQ - LIO + KDF + KXF -LPQ+ KXA + EPH + EXC + 

2LC - 2KF - 2ET - LDM - LXC + K1U + KXJ - LXT - LPO - LXM - 
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BIR - BPR - LIN - BIS + EPN - BXA - LPH + KPV + KDA - LDT - 
LXG - BXJ + KIS - BDJ + EXG + EXT - 2KJ + 2LT + KDW- BXF - 
BPV- BIU - BPU - BIV + KXW + 2BV+2BU + 2BR - 2EN + 2LH + 



KPR + EDT- BDF - 2KS - 2EC - 2KR + EIQ + EPQ + EPO + EIN + 
EDM - LDC + KDJ + EIH + EDG + 2BF + EDC + 2BA ) x 






(- 2xua^ + + lumQO^b + Izm^a^b + 



IxmjQ^b + 3mja^b^ + lutTijO^b + lumfib^ + lurriQab^ ) J 



The variables A ...X in equation (52) are symbols which are represented as 



A = e 



— {luniQb + lum^a + /7i^a6 + + jcm^a + + um^b) 



B = er/ 



fu + nija\ 
aj2 j 



C = e 



*■2 

{Aunifb + Axrrija + 3mQab + lumQO + IxmQO + + luniQb) 



D = 

£ = erf 



X — u + niQb 



bj2 



r-u + mQa'' 



aj2 



(52) 

(53) 

(54) 

(55) 

(56) 

(57) 
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(58) 



F = e 



ab 



(lurriQb + 2xm,a + ttigab + uniQa + xm^a + m]ab + um,b) 



G = e 



+ lumga + Ixm^a + m^,ab + lum^b) 



(59) 



H = e 



1 2 2 
— (2umjb + 2xmQa + rngOb^ nifax+ unifa + nijab + uniQb) 



1 = erf\ 



'X- u + mjb^ 



J = e 



1 2 2 
— ( 2umfO + + u/n^a + m^ojc + //i; 



= erf 



'-U + mja\ 



ajl 



(60) 

(61) 

(62) 

(63) 



L = erf 



fu + mQa^ 



aj2 



J 



(64) 



M = e 



*^2 2 

+ 4w/n^ + ^rrigab + 2uniQa + 2x/n^a + aZ? + 2umQb) 



(65) 



TV = e 



— (2umjb + 2uniQa + /n^a^ + /n;Ox + a/? + 



(66) 



2 2 

— (2w/ngfl + /ngai> + /n,ar + um,a + m,ab + um^b) 

O - e'*® 

fz-u- m,b\ 

' ■ ’"(-ir) 



(67) 

( 68 ) 



Q = e 



— (^2xmQa + mQab + m,ax+um,a+m^,ab + umQb) 



(69) 



/? = e 



^2 2 

(4jc/ngfl + mQab + 2xmfl + 3/n^ + 2unija + 2ump) 



(70) 
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(71) 



S = e 



1 2 2 

{Aum^a + nifyab + 2xm,a + 3m. ab + 2umja + 2umjb) 
2ab ^ ^ 



A 2 2 

-—7 {Aumja^3mQab^2umQa-¥2xmQa^mjab^2umQb) 



r = e 



t/ = e 



X 2 2 

{AuniQb + AxniQa + rngOft + 2jrm,a + 3miab + 2umfi + lumjb) 



V = e 



1 2 2 

(^umQb + + 3/n, + 2wm;<3 + 2umjb) 



1 



W = e 



— ( 2 xm /2 + mQab+ umQa-hmQOX + m^ab + um^b) 



X = erf 



fx — u — niQb^ 



bj2 



(72) 

(73) 

(74) 

(75) 

(76) 



a. Error Function Approximation 

The density function of the envelope detection approximations:^ given by 
(52) cannot be obtained in a mathematical closed form. The specific difficulty in doing this 
is to integrate expressions of the form 

^{erf{au-\-b)e‘^'^^*‘^'**^}du ( 77 ) 



or 

J {er/(au + b) erf{cu + d) du . 

A few closed form solutions exist for indefinite integrals of an error function multiplying 
an exponential, but their variable arguments are simpler and different in form from those 

of (77). Very accurate series expansions offering error residuals to the order of 10’^ to 



35 



10'^® also exist for the error function. These expressions, however, are of fifth or higher 
order, and are thus very difficult, if not impossible, to integrate in a closed form. [Ref. 9] 
The simplest approach available to integrate (52) in closed form is to 
approximate the error function with a polynomial. A second order polynomial-fitted 
approximation will ensure a closed form solution for the integration of terms in the form 
of (77) and (78). The error function can thus be approximated with a maximum residual 
error magnitude of less than 1.75e-3. In order to achieve this order of residual, the error 
function erfipc) has to be modelled in six segments. For the first five segments the error 
function is approximated by polynomials of the form 

erf{x) + + c (79) 

where the coefficients a, b and c for each segment are outlined in Table 2. For the sixth 
TABLE 2: ERROR FUNCTION APPROXIMATION POLYNOMIAL COEFFICIENTS 



segment 


a 


b 


c 


1: 0.0-0.5 


-0.2539 


1.1755 


-0.0018 


2; 0.5- 1.0 


-0.4739 


1.3547 


-0.0384 


3: 1.0- 1.5 


-0.2959 


0.9828 


0.1568 


4: 1. 5-2.0 


-0.0950 


0.3888 


0.5972 


5; 2.0-3.0 


-0.0075 


0.0409 


0.9440 


6: 3.0-®° 


0.0 


0.0 


1.0 



segment, the value of the error fimction is set to unity, since by constraint the error function 
is very close to a value of one for an argument value beyond 3.0. The comparison of the 
actual error function to its approximations, along with the corresponding residuals, are 
presented in Figures 5 to 14. Note that smaller residuals can be obtained by approximating 
with a higher order of polynomial. A conditional statement is included inside the 
integration loop to verify if the error function argument is negative. This is in order to limit 
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magnitude of residual 




Figure 5: Segment 1 of Error Function Approximation 




Figure 6: Residual of Segment 1 of Error Function Approximation 
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Figure 7: Segment 2 of Bror Function Approximation 




Figure 8: Residual of Segment 2 of Error Function Approximation 
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magnitude of residual 



0.98 




Figure 9: Segment 3 of Error Function Approximation 




Figure 10: Residual of Segment 3 of Error Function Approximation 
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magnitude ol residual 




Figure 11: Segment 4 of Error Function Approximation 




Figure 12: Residual of Segment 4 of Error Function Approximation 
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magnitude of residual 




Figure 13: Segment 5 of Error Fimction Approximation 




Figure 14: Residual of Segment 5 of Error Function Approximation 
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the number of additional integration terms introduced by the error function approximation. 
If the argument is negative, the value of the error function is found from the relation 
erf(—x) =—erf(x) [Ref. 9]. These curves are generated using the program in Appendix 
A. 

As an example of how to integrate (52) using the error function 
approximations, the first term (-2EQ) will be integrated here, for the range where the 
argument of the error function of E (57) is between 0.5 and 1.0 (segment 2). In other words, 
when the integration of (52) (integrating over u) is being performed between the limits 

u = a {rriQ - 1/J2) and u = a (niQ - J2 ) , integration of the first term (referred to as 
7i(z)) is expressed as 

a (ffig - J2) — 1 22ud‘ + + lum^a^b + 

7,(z) = J -2EQe^^^^ 

a(mQ- IZ-Jl) 

IzniQa^b + 2zm,ab + 3nija^b^ + 2um,a^b + 2um^ab^ + 2umQob ^ ) 



/ /o\ , Z~u + m„a\ 

a(ffig-V2) r _ 0.4739 [ -^1 + 1.3547 [ -0.0384 

_ _2 j J . { aj2 J { aj2 ) 

a(mQ- 1/72) 



^2 2 

— (IzniQa + niQab + znija + um^a + 



X 



— ^ (- Izua^ + + 3m^a^b^ + z^a^ + + 2um^a^b + Izmia^b^ 

^la^b^ C Q Q 



Iznijab + + lumjQ^b + 2umjab^ + 2umQab^ ) 



}du 
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or 



/l(z) = 



± 0 2 

0 00005 If ( 22 W*g + 2/ngi) + zm,+m;amg- lAl42m,a + m,b- 1.4142/ngfc) 
^ 

{mia + niQb)^ 



ba ( 16848mgZ?^ + 16848w;a^ + 5754.435a6m; - 9478Z>^ + 33696niiamQb + 

5754.435 mg ) ] - 



0.000025 ^2b 

(nifa + niQb) 



1 2 2 
— (42m^ + 4/n^Z> + 2zm/ + 2m/a/72g- L4142m;fl+ 2w;/?- 1.4142m^/?) 



6a (20819/ng6^ + 20819m^a^ + 24912.786a6/n^ - 189566^ + 41638m/a/ng6 + 



24912.7866^mQ ) ] 



(80) 



Using this error function approximation in (52) would result in a closed 
form target test cell density function (x) being composed of over 1 100 exponential 

terms. Hence, this may be a very costly method of determining (x) for a given situation. 

B. CLOSED FORM PDF BY CURVE-FITTING 

The probability density function obtained in Section ni.A.6 for the envelope 
detection approximation = aMax {\I\,\Q\} + bMin {\I\,\Q\} cannot be practically 
implemented for either system performance analysis or actual operational use. This is due 
to the very complex and tedious procedure required to convert the convolution integral of 
the test cell (52) into a closed form expression. Of course, the integral could be evaluated 
numerically, but this approach would be extremely computationally intensive and time- 
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consuming. This is especially true for the computation of probability of detection curves 
for the GO CFAR processor, which require at least four integrations (see Section III.A). 

The alternative presented here provides a relatively simple closed form expression for 
the density function of the target test cell with arbitrary SNR. Generally, the closed form 
solution for a target test cell density function can be modelled as an exponential of the form 



where the coefficients A, B, ..... K are functions of the SNR and may be expressed as 
polynomials of the form 

A - Ag/ + A’jS^ + Ag/ +A^s^+ A 4 / + A'^s^ + ^ 2 ^^ + A + Aq (82) 

with similar expansions for B , ..... K. The variable s represents the SNR in decibels. 

Polynomial curve-fitting in a least-squares sense is used to yield (81). Such curve- 
fitting, using a tenth order exponential as shown above, generally provides a maximum 
residual error magnitude of 0.03 if the range of the density function over which the curve- 
fitting is applied is small enough. Now, as is well known, the constraint of “small enough” 
is ambiguous. In this research, the range over which the density function is greater than 

10'^ was considered the limit of the curve-fitting range. A tenth order exponential is used 
in (8 1) so as to maintain a minimiun residual error. These, however, have the disadvantage 
of not being integrable in closed form, and have to be integrated numerically. Higher order 
exponentials provided no fiuther improvement in the residual error. Smaller order 
exponentials in (81) would ensure that a closed form solution could be obtained for the 
integration of the target test cell density function p<-(x) in (7). Such smaller order 
exponentials, however, proved to be very inaccurate. 
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The solution of the form (81) is obtained by performing the numerical convolution 
(50) for a number of increasing signal-to-noise ratios and for a given combination of the 
multiplying constants a and b of the envelope detection approximation. Each numerical 
convolution was averaged over the phase angle in the range 0 ...k/ 4, as described in 
Section n.B. The results of such an operation are shown in Figure 15, where the 



0.5r 




Figure 15: Test Cell PDF for Increasing SNR 



consecutive density function curves correspond to the test cell density for each of the 
signal-to-noise ratios (SNRs) used and where the constants were chosen to be a=1.0 and 
Z?=0.25. These particular constant values were chosen arbitrarily in order to demonstrate the 
present discussion. These curves are generated using the program in Appendix B. 
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Each curve of Figure 15 is then curve-fitted, independently of the others, in the form 
of (81). For example, modeling of the first curve (SNR = 0 dB) over the range 

X = 0.01 ...5.6 (region where > 10“^) yields the curve-fit shown in Figure 16, which 



06 



05 



04 
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0.2 



0.1 




3 

X 



Figure 16: Curve- Fit of (x) for SNR = 0 dB, a = 1.0, 6 = 0.25 



also shows the corresponding magnitude of the residual. The exponential coefficients 
required for this curve-fit are 



A = -0.0015 



B = 0.0459 
C = -0.5889 



D = 4.262 
E = -19.1293 
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F = 55.1303 



G = -102.2169 
H = 119.0197 
7 = -83.2403 
J = 33.1493 
a: = -7.482 



Outside of the curve-fit range, the value of {x) is taken to be zero. 

By performing a number of these curve-fits for successive (x) curves in Figure 15, 

sequences of coefficients for each of K can be plotted out. AH the curve-fits thus 

obtained are shown in Figure 17, with their corresponding residual error magnitudes. The 
coefficient sequences are plotted out for SNR = 0....21 dB and for a = 1.0 and b = 0.25 in 
Figures 18 to 28. It is rather difficult to extract valid coefficient values from these plots. 
Therefore the exact coefficient values used to produce the plots are tabulated. This is 
actually done for all the envelope detection approximation situations outlined in Table 1 so 
as to provide a good reference design tool. One density p^ (x) can be produced at each 
integer value of SNR. Tables 3 to 9 provide all the necessary coefficients. The SNR range 
0....21 dB was chosen since it covers the range over which Pj) is less than one (see section 
V). Beyond this range, Pj) is constant at one. 

The curves for A^, ,K can in turn be curve-fitted in the form of (82). For instance, 

the exponential coefficient A of (81) can be expressed as (82) where the coefficients 

Ay Aq are determined by the polynomial curve-fit. Thus, for a given a and b and SNR 

range, it appears that the density function of a target test cell could be approximated by both 
(81) and (82). This would, in effect, remove the dependence of (81) on the SNR. The 
coefficient curves are also compared to their curve-fits in Figures 18 to 28. As well, these 
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0.5^ A- SNR = 0,1 ,21dB 
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X 



Figure 17: Curve-Fits for (jt) , with Residual Errors 
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Figure 18: Behaviour of Exponential Coefficient A as a Function of SNR 




Figure 19: Behaviour of Exponential Coefficient B as a Function of SNR 
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Figure 20: Behaviour of Exponential Coefficient C as a Function of SNR 




Figure 21: Behaviour of Exponential Coefficient D as a Function of SNR 
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Figure 22: Behaviour of Exponential Coefficient E as a Function of SNR 




Figure 23: Behaviour of Exponential Coefficient F as a Function of SNR 



51 




Figure 24: Behaviour of Exponential Coefficient G as a Function of SNR 




Figure 25: Behaviour of Exponential Coefficient H as a Function of SNR 



52 




SNR. dB 



Figure 26: Behaviour of Exponential Coefficient I as a Function of SNR 




SNR. dB 



Figure 27: Behaviour of Exponential Coefficient J as a Function of SNR 
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Figure 28: Behaviour of Exponential Coefficient K as a Function of SNR 
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show the corresponding residual error magnitudes. Note that the residual error in each case 
hovers very near zero, indicating an effective curve-fit. As it turns out, however, attempts 
to obtain a reasonable target test cell density function using the curve-fits of the coefficients 

A,B K were unsuccessful. It appears that the result of a polynomial expansion as in (82) 

is very sensitive to slight changes in its coefficient values. Hence, since a curve-fit 
generally has errors at every point, regardless of how small, the resultant polynomial does 
not behave as expected. 

As mentioned earlier, smaller order exponential expansions approximating {x) 
generate significant errors, as can be seen in Figure 29, which shows a third order 
exponential curve-fit. 

The program used to generate all the curve-fits and plots in this discussion is in 
Appendix C. 
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TABLE 3: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 1.0, FOR 



SNR (dB) 


Range 
Lower 
Limit Xl 


Range 
Upper 
Limit Xy 


A 


B 


C 


0 


.02 


6.78 


-1.9238e-4 


6.9423e-3 


-1.0802e-l 


1 


.025 


7.02 


-1.2619e-4 


4.7207e-3 


-7.6157e-2 


2 


.03 


7.295 


-8.0869e-5 


3.1468e-3 


-5.2814e-2 


3 


.04 


7.6 


-4.7671e-5 


1.9360e-3 


-3.3925e-2 


4 


.055 


7.945 


-2.6193e-5 


1.1146e-3 


-2.0473e-2 


5 


.08 


8.33 


-1.3038e-5 


5.8365e-4 


-1.1284e-2 


6 


.125 


8.765 


-5.5545e-6 


2.6303e-4 


-5.3839e-3 


7 


.22 


9.25 


-1.7716e-6 


8.9455e-5 


-1.9557e-3 


8 


.4 


9.8 


-3.9263e-7 


2.1412e-5 


-5.0709e-4 


9 


.705 


10.415 


-5.3040e-8 


3.2497e-6 


-8.7064e-5 


10 


1.11 


11.105 


2.3166e-8 


-1.2966e-6 


3.0154e-5 


11 


1.61 


11.88 


4.6389e-8 


-2.9893e-6 


8.2487e-5 


12 


2.2 


12.745 


4.9775e-8 


-3.5120e-6 


1.0648e-4 


13 


2.89 


13.72 


4.5677e-8 


-3.5275e-6 


1.1728e-4 


14 


3.68 


14.81 


4.0119e-8 


-3.4077e-6 


1.2490e-4 


15 


4.58 


16.035 


3.1905e-8 


-2.9665e-6 


1.1892e-4 


16 


5.59 


17.405 


2.1449e-8 


-2.1390e-6 


9.0897e-5 


17 


6.715 


18.945 


1.0642e-8 


-1.0544e-6 


4.1129e-5 


18 


7.97 


20.67 


1.8897e-9 


3.6074e-9 


-1.8714e-5 


19 


9.375 


22.605 


-3.8683e-9 


8.4043e-7 


-7.6459e-5 


20 


10.945 


24.775 


-6.3980e-9 


1.3038e-6 


-1.1726e-4 


21 


12.7 


27.21 


-6.5285e-9 


1.4038e-6 


-1.3443e-4 
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TABLE 3 (cont’d): EXPONENTIAL COEFFICIENTS, o = 1,0. ^ = 1.0, FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


2.000e-2 


6.7800e+0 


-1.9238e-4 


6.9423e-3 


-1.0802e-l 


1 


2.500e-2 


7.0200e+0 


-1.2619e-4 


4.7207e-3 


-7.6157e-2 


2 


3.000e-2 


7.2950e+0 


-8.0869e-5 


3.1468e-3 


-5.2814e-2 


3 


4.000e-2 


7.6000e+0 


-4.767 le-5 


1.9360e-3 


-3.3925e-2 


4 


5.500e-2 


7.9450e+0 


-2.6193e-5 


1.1146e-3 


-2.0473e-2 


5 


8.000e-2 


8.3300e+0 


-1.3038e-5 


5.8365e-4 


-1.1284e-2 


6 


1.250e-l 


8.7650e+0 


-5.5545e-6 


2.6303e-4 


-5.3839e-3 


7 


2.200e-l 


9.2500e+0 


-1.7716e-6 


8.9455e-5 


-1.9557e-3 


8 


4.000e-l 


9.8000e+0 


-3.9263e-7 


2.1412e-5 


-5.0709e-4 


9 


7.050e-l 


1.0415e+l 


-5.3040e-8 


3.2497e-6 


-8.7064e-5 


10 


l.llOOe+0 


1.1105e+l 


2.3166e-8 


-1.2966e-6 


3.0154e-5 


11 


1.6100e+0 


1.1880e+l 


4.6389e-8 


-2.9893e-6 


8.2487e-5 


12 


2.2000e+0 


1.2745e+l 


4.9775e-8 


-3.5120e-6 


1.0648e-4 


13 


2.8900e+0 


1.3720e+l 


4.5677e-8 


-3.5275e-6 


1.1728e-4 


14 


3.6800e+0 


1.4810e+l 


4.0119e-8 


-3.4077e-6 


1.2490e-4 


15 


4.5800e+0 


1.6035e+l 


3.1905e-8 


-2.9665e-6 


1.1892e-4 


16 


5.5900e+0 


1.7405e+l 


2.1449e-8 


-2.1390e-6 


9.0897e-5 


17 


6.7150e+0 


1.8945e+l 


1.0642e-8 


-1.0544e-6 


4.1129e-5 


17 


7.9700e+0 


2.0670e+l 


1.8897e-9 


3.6074e-9 


-1.8714e-5 


19 


9.3750e+0 


2.2605e+l 


-3.8683e-9 


8.4043e-7 


-7.6459e-5 


20 


1.0945e+l 


2.4775e+l 


-6.3980e-9 


1.3038e-6 


-1.1726e-4 


21 


1.2700e+l 


2.7210e+l 


-6.5285e-9 


1.4038e-6 


-1.3443e-4 
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TABLE 3 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, 6 = 1.0, FOR ;c^ 



SNR (dB) 


I 


J 


K 


0 


-5.000e+l 


2.568 le+1 


-8.277e+0 


1 


-4.434e+l 


2.4021e+l 


-8.463e+0 


2 


-3.929e+l 


2.2521e+l 


-8.740e+0 


3 


-3.330e+l 


2.0556e+l 


-9.04 le+0 


4 


-2.733e+l 


1.8463e+l 


-9.418e+0 


5 


-2.129e+l 


1.6166e+l 


-9.875e+0 


6 


-1.528e+l 


1.3658e+l 


-1.042e+l 


7 


-9.420e+0 


1.0912e+l 


-1.106e+l 


8 


-5.048e+0 


. 8.6388e+0 


-1.196e+l 


9 


-2.679e+0 


7.3700e+0 


-1.332e+l 


10 


-1.359e+0 


6.7302e+0 


-1.520e+l 


11 


-1.7187e-2 


5.9774e+0 


-1.748e+l 


12 


1.4997e+0 


5.1285e+0 


-2.041e+l 


13 


3.415 le+0 


4.4415e+0 


-2.482e+l 


14 


6.7708e+0 


3.875 le+0 


-3.266e+l 


15 


3.2824e+0 


2.4171e+l 


-6.785e+l 


16 


-4.670e+l 


1.7287e+2 


-2.608e+2 


17 


-2.556e+2 


8.0545e+2 


-1.105e+3 


18 


-8.599e+2 


2.8178e+3 


-4.058e+3 


19 


-2.293e+3 


8.1620e+3 


-1.284e+4 


20 


-5.103e+3 


2.0022e+4 


-3.484e+4 


21 


-9.761e+3 


4.2408e+4 


-8.194e+4 
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TABLE 4: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5, FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xjj 


A 


B 


C 


0 


.015 


5.89 


-8.2002e-4 


2.5679e-2 


-3.4668e-l 


1 


.015 


6.085 


-5.9881e-4 


1.9372e-2 


-2.7018e-l 


2 


.02 


6.3 


-3.8511e-4 


1.2918e-2 


-1.8686e-l 


3 


.03 


6.545 


-2.2335e-4 


7.8035e-3 


-1.1764e-l 


4 


.04 


6.82 


-1.3015e-4 


4.7479e-3 


-7.4765e-2 


5 


.055 


7.125 


-7.1127e-5 


2.7179e-3 


-4.4848e-2 


6 


.09 


7.47 


-3.1612e-5 


1.2734e-3 


-2.2174e-2 


7 


.155 


7.86 


-1.1488e-5 


4.9103e-4 


-9.0876e-3 


8 


.3 


8.29 


-2.763 le-6 


1.2652e-4 


-2.5165e-3 


9 


.57 


8.775 


-4.8343e-7 


2.3765e-5 


-5.0945e-4 


10 


.955 


9.32 


-1.1659e-7 


5.9043e-6 


-1.2978e-4 


11 


1.44 


9.93 


-3.4433e-8 


1.6178e-6 


-3.0537e-5 


12 


2.02 


10.615 


-1.1632e-8 


3.8520e-7 


1.3429e-7 


13 


2.69 


11.38 


-1.0728e-8 


4.5402e-7 


-3.8989e-6 


14 


3.465 


12.235 


-1.3047e-8 


7.5926e-7 


-1.6489e-5 


15 


4.35 


13.2 


-1.0466e-8 


6.6902e-7 


-1.5785e-5 


16 


5.345 


14.275 


-5.4726e-9 


3.0051e-7 


-2.3702e-6 


17 


6.47 


15.485 


9.112e-10 


-3.1705e-7 


2.6409e-5 


18 


7.73 


16.845 


7.4178e-9 


-1.0975e-6 


7.0533e-5 


19 


9.145 


18.365 


1.2680e-8 


-1.8827e-6 


1.2429e-4 


20 


10.720 


20.075 


1.4164e-8 


-2.246 le-6 


1.5917e-4 


21 


12.485 


21.99 


9.0439e-9 


-1.5160e-6 


1.1354e-4 
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TABLE 4 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, ^) = 0.5, FOR jc^ 



SNR(dB) 


D 


E 


F 


G 


H 


0 


2.6410e+0 


-1.248e+l 


3.7896e+l 


-7.414e+l 


9.1405e+l 


I 


2.1262e+0 


-1.038e+l 


3.2549e+l 


-6.575e+l 


8.3637e+l 


2 


1.5258e+0 


-7.732e+0 


2.5194e+l 


-5.293e+l 


7.0114e+l 


3 


1.0018e+0 


-5.299e+0 


1.8047e+l 


-3.970e+l 


5.521 le+1 


4 


6.6539e-l 


-3.68 le+0 


1.3122e+l 


-3.025e+l 


4.4157e+l 


5 


4.1852e-l 


-2.430e+0 


9.0994e+0 


-2.207e+l 


3.3975e+l 


6 


2.1867e-l 


-1.344e+0 


5.3428e+0 


-1.381e+l 


2.277 le+1 


7 


9.545 le-2 


-6.2675e-l 


2.6724e+0 


-7.451e+0 


1.3385e+l 


8 


2.8529e-2 


-2.0342e-l 


9.501 le-1 


-2.939e+0 


5.9730e+0 


9 


6.2662e-3 


-4.9004e-2 


2.5533e-l 


-9.0354e-l 


2.1830e+0 


10 


1.6361e-3 


-1.3242e-2 


7.343 le-2 


-2.9130e-l 


8.5036e-l 


11 


2.875 le-4 


-1.3054e-3 


1.7862e-3 


5.5201e-4 


6.7262e-2 


12 


-1.8672e-4 


3.7803e-3 


-3.6964e-2 


2.0507e-l 


-6.4428e-l 


13 


-1.2208e-4 


3.6073e-3 


-4.3359e-2 


2.9147e-l 


-1.150e+0 


14 


1.2564e-4 


1.0598e-3 


-3.1541e-2 


3.0413e-l 


-1.580e+0 


15 


1.1806e-4 


1.8113e-3 


-5.1991e-2 


5.6644e-l 


-3.452e+0 


16 


-2.2027e-4 


8.3078e-3 


-1.4429e-l 


1.4839e+0 


-9.526e+0 


17 


-1.0750e-3 


2.5944e-2 


-4.0165e-l 


4.1088e+0 


-2.769e+l 


18 


-2.6083e-3 


6.1696e-2 


-9.7786e-l 


1.0538e+l 


-7.635e+l 


19 


-4.805 le-3 


1.2050e-l 


-2.049e+0 


2.3915e+l 


-1.893e+2 


20 


-6.6375e-3 


1.8039e-l 


-3.339e+0 


4.2635e+l 


-3.709e+2 


21 


-5.0050e-3 


1.4385e-l 


-2.818e+0 


3.8142e+l 


-3.524e+2 
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TABLE 4 (cont’d): EXPONENTIAL COEFFICIENTS, a=\.0,b = 0.5, FOR ;c^ 



SNR (dB) 


I 


J 


K 


0 


-6.837e+l 


3.0021e+l 


-7.912e+0 


1 


-6.447e+l 


2.928 le+1 


-8.219e+0 


2 


-5.646e+l 


2.7182e+l 


-8.458e+0 


3 


-4.696e+l 


2.4447e+l 


-8.708e+0 


4 


-3.961e+l 


2.2285e+l 


-9.119e+0 


5 


-3.235e+l 


1.9990e+l 


-9.63 le+0 


6 


-2.360e+l 


1.6866e+l 


-1.016e+l 


7 


-1.553e+l 


1.3674e+l 


-1.083e+l 


8 


-8.343e+0 


1.0434e+l 


-1.164e+l 


9 


-4.097e+0 


8.3157e+0 


-1.289e+l 


10 


-2.318e+0 


7.4870e+0 


-1.476e+l 


11 


-1.013e+0 


6.8070e+0 


-1.709e+l 


12 


5.1247e-l 


5.5877e+0 


-1.964e+l 


13 


2.0785e+0 


3.7456e+0 


-2.230e+l 


14 


4.1725e+0 


-2.4275e-2 


-2.386e+l 


15 


1.1826e+l 


-1.600e+l 


-1.509e+l 


16 


3.7092e+l 


-7.459e+l 


3.7069e+l 


17 


1.1787e+2 


-2.81 le+2 


2.6044e+2 


18 


3.5586e+2 


-9.573e+2 


1.0984e+3 


19 


9.725 le+2 


-2.919e+3 


3.8492e+3 


20 


2.1027e+3 


-7.009e+3 


1.0380e+4 


21 


2.1290e+3 


-7.587e+3 


1.2057e+4 
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TABLE 5; EXPONENTIAL COEFFICIENTS, a=l.O,b = 0.25, FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.01 


5.6 


-1.5422e-3 


4.5899e-2 


-5.8885e-l 


1 


.01 


5.77 


-1.1524e-3 


3.5342e-2 


-4.6724e-l 


2 


.015 


5.965 


-7.3691e-4 


2.3410e-2 


-3.2073e-l 


3 


.02 


6.185 


-4.6829e-4 


1.545 le-2 


-2.1993e-l 


4 


.025 


6.435 


-2.9259e-4 


1.0057e-2 


-1.4916e-l 


5 


.04 


6.715 


-1.5489e-4 


5.5750e-3 


-8.6645e-2 


6 


.065 


7.025 


-7.4596e-5 


2.8221e-3 


-4.6143e-2 


7 


.11 


7.375 


-3.0750e-5 


1.2297e-3 


-2.1284e-2 


8 


.22 


7.77 


-8.6030e-6 


3.6777e-4 


-6.8269e-3 


9 


.45 


8.21 


-1.3960e-6 


6.4832e-5 


-1.3188e-3 


10 


.815 


8.705 


6.4795e-8 


-2.7352e-6 


4.0778e-5 


11 


1.285 


9.255 


2.3171e-7 


-1.2180e-5 


2.7321e-4 


12 


1.84 


9.88 


2.4897e-7 


-1.4384e-5 


3.5862e-4 


13 


2.48 


10.575 


2.2852e-7 


-1.4434e-5 


3.9556e-4 


14 


3.205 


11.360 


1.7959e-7 


-1.2307e-5 


3.6627e-4 


15 


4.03 


12.240 


1.0541e-7 


-7.5467e-6 


2.307 le-4 


16 


4.955 


13.225 


7.761 7e-9 


5.1223e-7 


-7.3435e-5 


17 


5.99 


14.335 


-1.0992e-7 


1.2479e-5 


-6.2665e-4 


18 


7.14 


15.58 


-2.3077e-7 


2.7493e-5 


-1.4599e-3 


19 


8.42 


16.975 


-3.3934e-7 


4.4222e-5 


-2.5740e-3 


20 


9.84 


18.545 


-4.1982e-7 


6.0504e-5 


-3.8995e-3 


21 


11.425 


20.310 


-4.6637e-7 


7.4716e-5 


-5.3582e-3 
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TABLE 5 (cont’d): EXPONENTIAL COEFFICIENTS, a=l.O,b = 0.25, FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


4.2620e+0 


-1.913e+l 


5.5130e+l 


-1.022e+2 


1.1902e+2 


1 


3.485 le+0 


-1.612e+l 


4.7874e+l 


-9.146e+l 


1.0967e+2 


2 


2.4803e+0 


-1.190e+l 


3.6714e+l 


-7.293e+l 


9.1111e+l 


3 


1.7678e+0 


-8.822c+0 


2.8317e+l 


-5.859c+l 


7.6320e+l 


4 


1.2496e+0 


-6.503e+0 


2.1771e+l 


-4.701e+l 


6.395 le+1 


5 


7.6135e-l 


-4.160e+0 


1.4649e+l 


-3.334e+l 


4.7930e+l 


6 


4.2708e-l 


-2.462e+0 


9.1652e+0 


-2.2llle+l 


3.3833e+l 


7 


2.0893e-l 


-1.28 le+0 


5.0835e+0 


-1.313e+l 


2.1641e+l 


8 


7.2173e-2 


-4.7894e-l 


2.0726e+0 


-5.884e+0 


1.0776e+l 


9 


1.5444e-2 


-1.1498e-l 


5.6590e-l 


-1.852e+0 


3.9729e+0 


10 


-1.4745e-4 


-2.8661e-3 


4.2063e-2 


-2.5654e-l 


8.7002e-l 


11 


-3.3979e-3 


2.5470e-2 


-1.1715e-l 


3.2201e-l 


-4.7101e-l 


12 


-5.043 le-3 


4.3919e-2 


-2.4549e-l 


8.8803e-l 


-2.052e+0 


13 


-6.1599e-3 


6.0012e-2 


-3.8036e-l 


1.5864e+0 


-4.316e+0 


14 


-6.1978e-3 


6.561 le-2 


-4.5101e-l 


2.0280e+0 


-5.873e+0 


15 


-3.8945e-3 


3.8882e-2 


-2.2216e-l 


5.5003e-l 


1.0567e+0 


16 


3.0544e-3 


-6.7722e-2 


9.2140e-l 


-8.070e+0 


4.5690e+l 


17 


1.8324e-2 


-3.4534e-l 


4.3796e+0 


-3.783e+l 


2.1960e+2 


18 


4.547 le-2 


-9.1937e-l 


1.2599e+l 


-1.185e+2 


7.5399e+2 


19 


8.808 le-2 


-1.962e+0 


2.9692e+l 


-3.093e+2 


2.1882e+3 


20 


1.4796e-l 


-3.659e+0 


6.1622e+l 


-7.153e+2 


5.6515e+3 


21 


2.2647e-l 


-6.246e+0 


1.1743e+2 


-1.524e+3 


1.3479e+4 
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TABLE 5 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, 6 = 0.25, FOR 



SNR (dB) 


I 


J 


K 


0 


-8.324e+l 


3.3149e+l 


-7.482e+0 


1 


-7.892e+l 


3.2466e+l 


-7.792e+0 


2 


-6.856e+l 


2.9934e+l 


-7.998e+0 


3 


-6.005e+l 


2.785 le+1 


-8.332e+0 


4 


-5.266e+l 


2.6023e+l 


■8.80c+0 


5 


-4.204e+l 


2.2912e+l 


-9.243e+0 


6 


-3.193e+l 


1.9693e+l 


-9.808e+0 


7 


-2.239e+l 


1.6344e+l 


-1.052e+l 


8 


-1.288e+l 


1.2518e+l 


-1.132e+l 


9 


-6.047e+0 


9.3673e+0 


-1.243e+l 


10 


-2.370e+0 


7.5056e+0 


-1.409e+l 


11 


-4.4474e-l 


6.5170e+0 


-1.637e+l 


12 


2.2908e+0 


4.4562e-fO 


-1.873e+l 


13 


6.8048e+0 


6.2160e-2 


-2.059e+l 


14 


9.8958e+0 


-2.284e+0 


-2.487e+l 


15 


-1.177e+l 


3.8159e+l 


-6.326e+l 


16 


-1.622e+2 


3.3524e+2 


-3.290e+2 


17 


-8.207e+2 


1.7913e+3 


-1.763e+3 


18 


-3.110e+3 


7.5127e+3 


-8.109e+3 


19 


-1.006e+4 


2.7171e+4 


-3.276e+4 


20 


-2.908e+4 


8.7994e+4 


-1.190e+5 


21 


-7.775e+4 


2.6410e+5 


-4.013e+5 
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TABLE 6: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8, FOR ;c^ 



SNR (dB) 


Range 
Lower 
Limit Xl 


Range 
Upper 
Limit Xu 


A 


B 


C 


0 


.015 


5.735 


-1.0551e-3 


3.2194e-2 


-4.2355e-l 


1 


.015 


5.915 


-7.8236e-4 


2.462 le-2 


-3.3409e-l 


2 


.02 


6.120 


-5.0488e-4 


1.6466e-2 


-2.3164e-l 


3 


.025 


6.350 


-3.2244e-4 


1.0926e-2 


-1.5975e-l 


4 


.035 


6.610 


-1.8673e-4 


6.6024e-3 


-1.0077e-l 


5 


.05 


6.900 


-1.0127e-4 


3.7489e-3 


-5.9940e-2 


6 


.08 


7.230 


-4.6655e-5 


1.8189e-3 


-3.0660e-2 


7 


.135 


7.595 


-1.8049e-5 


7.4499e-4 


-1.3316e-2 


8 


.265 


8.005 


-4.5077e-6 


1.9928e-4 


-3.8292e-3 


9 


.515 


8.465 


-7.2031e-7 


3.4402e-5 


-7.1993e-4 


10 


.895 


8.980 


-6.0718e-8 


2.9103e-6 


-6.3560e-5 


11 


1.375 


9.560 


6.8766e-8 


-4.2046e-6 


1.0972e-4 


12 


1.945 


10.205 


1.0133e-7 


-6.5844e-6 


1.8570e-4 


13 


2.605 


10.93 


1.0839e-7 


-7.6542e-6 


2.3624e-4 


14 


3.365 


11.745 


1.0667e-7 


-8.2503e-6 


2.8031e-4 


15 


4.230 


12.660 


1.0218e-7 


-8.6974e-6 


3.2652e-4 


16 


5.210 


13.685 


9.3045e-8 


-8.7357e-6 


3.6289e-4 


17 


6.310 


14.835 


7.6490e-8 


-7.9157e-6 


3.6312e-4 


18 


7.555 


16.125 


5.1897e-8 


-5.8704e-6 


2.9418e-4 


19 


8.945 


17.575 


2.1102e-8 


-2.4552e-6 


1.2384e-4 


20 


10.505 


19.20 


-1.2266e-8 


2.1881e-6 


-1.7090e-4 


21 


12.255 


21.03 


-4.2599e-8 


7.4224e-6 


-5.7898e-4 
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TABLE 6 (cont’d): EXPONENTIAL COEFFICIENTS, o = 1.0, b = 3/8, FOR 



SNR(dB) 


D 


E 


F 


G 


H 


0 


3.1452e+0 


-1.449e+l 


4.2928e+l 


-8.197e+l 


9.8681c+l 


1 


2.5587e+0 


-1.216e+l 


3.7145e+l 


-7.312e+l 


9.0702e+l 


2 


1.8401e+0 


-9.075e+0 


2.8793e+l 


-5.894e+l 


7.6128e+l 


3 


1.3193e+0 


-6.768e+0 


2.2347e+l 


-4.764e+l 


6.4139e+l 


4 


8.6934e-l 


-4.662e+0 


1.6111e+l 


-3.60e+l 


5.0895e+l 


5 


5.4208e-l 


-3.05 le+0 


1.1077e+l 


-2.605e+l 


3.8865e+l 


6 


2.927 le-1 


-1.742e+0 


6.7053e+0 


-1.677e+l 


2.6737e+l 


7 


1.3509e-l 


-8.5685e-l 


3.5282e+0 


-9.489e+0 


1.6392e+l 


8 


4.1968e-2 


-2.8944e-l 


1.3072e+0 


-3.90e+0 


7.5929e+0 


9 


8.6913e-3 


-6.7027e-2 


3.448 le-1 


-1.197e+0 


2.7836e+0 


10 


8.5778e-4 


-8.0720e-3 


5.4730e-2 


-2.6132e-l 


8.4939e-l 


11 


-1.5920e-3 


1.4004e-2 


-7.6259e-2 


2.503 le-1 


-4.3658e-l 


12 


-2.9748e-3 


2.9754e-2 


-1.9258e-l 


8.1016e-l 


-2.164e+0 


13 


-4.1790e-3 


4.6704e-2 


-3.4300e-l 


1.6704e+0 


-5.314e+0 


14 


-5.4932e-3 


6.8557e-2 


-5.6782e-l 


3.1554e+0 


-1.161e+l 


15 


-7.1065e-3 


9.9104e-2 


-9.2385e-l 


5.8253e+0 


-2.454e+l 


16 


-8.7713e-3 


1.3643e-l 


-1.426e+0 


1.0129e+l 


-4.833e+l 


17 


-9.7131e-3 


1.6759e-l 


-1.948e+0 


1.5427e+l 


-8.226e+l 


18 


-8.5878e-3 


1.6147e-l 


-2.040e+0 


1.7502e+l 


-1.006e+2 


19 


-3.5170e-3 


6.0717e-2 


-6.2883e-l 


3.2831e+0 


1.5440e+0 


20 


7.7315e-3 


-2.2503e-l 


4.4116e+0 


-5.908e+l 


5.3408e+2 


21 


2.6618e-2 


-7.9853e-l 


1.6329e+l 


-2.304e+2 


2.2156e+3 
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TABLE 6 (cont’d): EXPONENTIAL COEFnCIENTS, a = 1.0, /> = 3/8, FOR 



SNR (dB) 


I 


J 


K 


0 


-7.202e+l 


3.0644e+l 


-7.6636-40 


1 


-6.813e+l 


2.9955e+l 


-7.9736-40 


2 


-5.974e+l 


2.7835e+l 


-8.2136-40 


3 


-5.259e+l 


2.600764-1 


-8.5706-40 


4 


-4.405e+l 


2.357964-1 


-8.9606-40 


5 


-3.575e+l 


2.104064-1 


-9.4536-40 


6 


-2.663e+l 


1.794364-1 


-1.00164-1 


7 


-1.809e+l 


1.473364-1 


-1.0706-41 


8 


-9.943e+0 


1.12266-41 


-1.15164-1 


9 


-4.812e+0 


8.75246-40 


-1.27164-1 


10 


-2.370e+0 


7.55146-40 


-1.45164-1 


11 


-3.7835e-l 


6.3716640 


-1.6686-41 


12 


2.9060e+0 


3.473964-0 


-1.85664-1 


13 


1.0025e+l 


-4.9036-40 


-1.7936-41 


14 


2.6601e+l 


-2.8816-41 


-7.431640 


15 


6.5662e+l 


-9.50764-1 


3.594964-1 


16 


1.4789e+2 


-2.5606-42 


1.6626642 


17 


2.8225e+2 


-5.5616-42 


4.486464-2 


18 


3.6907e+2 


-7.7146-42 


6.56516-42 


19 


-1.287e+2 


7.037164-2 


-1.3506-43 


20 


-3.123e+3 


1.0677644 


-1.627644 


21 


-1. 389644 


5.1234644 


-8.455644 
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TABLE 7: EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR jc^ 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.01 


5.575 


-1.5822e-3 


4.6842e-2 


-5.9773e-l 


1 


.015 


5.75 


-1.0277e-3 


3.1439e-2 


-4.1472e-l 


2 


.02 


5.95 


-6.6134e-4 


2.097 le-2 


-2.8683e-l 


3 


.025 


6.175 


-4.2092e-4 


1.3871e-2 


-1.9723e-l 


4 


.03 


6.43 


-2.6300e-4 


9.0345e-3 


-1.3393e-l 


5 


.045 


6.715 


-1.3964e-4 


5.0252e-3 


-7.8090e-2 


6 


.075 


7.03 


-6.3592e-5 


2.4088e-3 


-3.9441e-2 


7 


.13 


7.39 


-2.3968e-5 


9.6210e-4 


-1.6722e-2 


8 


.255 


7.79 


-6.0010e-6 


2.5797e-4 


-4.8188e-3 


9 


.5 


8.235 


-9.5220e-7 


4.4183e-5 


-8.9789e-4 


10 


.87 


8.74 


-9.327e-8 


4.3495e-6 


-9.1379e-5 


11 


1.33 


9.30 


7.3921e-8 


-4.5278e-6 


1.1780e-4 


12 


1.885 


9.93 


1.1898e-7 


-7.6182e-6 


2.1161e-4 


13 


2.525 


10.635 


1.2892e-7 


-8.9440e-6 


2.7120e-4 


14 


3.26 


11.43 


1.2802e-7 


-9.7162e-6 


3.2398e-4 


15 


4.10 


12.315 


1.2580e-7 


-1.0500e-5 


3.8669e-4 


16 


5.055 


13.315 


1.1921e-7 


-1.0978e-5 


4.4761e-4 


17 


6.125 


14.435 


1.0417e-7 


-1.0596e-5 


4.7828e-4 


18 


7.335 


15.69 


7.995 le-8 


-8.9591e-6 


4.4593e-4 


19 


8.69 


17.1 


4.7268e-8 


-5.7463e-6 


3.0943e-4 


20 


10.205 


18.685 


9.6952e-9 


-9.9316e-7 


3.6398e-5 


21 


11.905 


20.46 


-2.6865e-8 


4.7589e-6 


-3.7602e-4 
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TABLE 7 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


4.3026e+0 


-1.920e+l 


5.5036e-t-l 


-1.015e+2 


1.1772e+2 


1 


3.0878e+0 


-1.427e+l 


4.2366e+l 


-8.109e+l 


9.7818e+l 


2 


2.2153e+0 


-1.062e+l 


3.2776e+l 


-6.524e+l 


8.1971e+l 


3 


1.5841e+0 


-7.903e+0 


2.5382e+l 


-5.264e+l 


6.895 le+1 


4 


1.1217e+0 


-5.838e+0 


1.9563e+l 


-4.235e+l 


5.7924e+l 


5 


6.8618e-l 


-3.751e+0 


1.3220e+l 


-3.016e+l 


4.3596e+l 


6 


3.6570e-l 


-2.113e+0 


7.8944e+0 


-1.916e+l 


2.9609e+l 


7 


1.6495e-l 


-1.017e+0 


4.0704e+0 


-1.064e+l 


1.7857e+l 


8 


5.1331e-2 


-3.4399e-l 


1.5092e+0 


-4.373e+0 


8.2716e+0 


9 


1.0521e-2 


-7.8728e-2 


3.9300e-l 


-1.325e+0 


2.9972e+0 


10 


1.1689e-3 


-1.0313e-2 


6.5609e-2 


-2.9716e-l 


9.2877e-l 


11 


-1.6974e-3 


1.4775e-2 


-7.9292e-2 


2.5466e-l 


-4.2467e-l 


12 


-3.337 le-3 


3.2841e-2 


-2.0899e-l 


8.6341e-l 


-2.259e+0 


13 


-4.7126e-3 


5.1728e-2 


-3.7300e-l 


1.7824e+0 


-5.557e+0 


14 


-6.2315e-3 


7.6337e-2 


-6.2056e-l 


3.3838e+0 


-1.221e+l 


15 


-8.2585e-3 


1.1305e-l 


-1.035e+0 


6.4082e+0 


-2.65 le+1 


16 


-1.0625e-2 


1.6243e-l 


-1.669e+0 


1.1671e+l 


-5.484e+l 


17 


-1.2603e-2 


2.1449e-l 


-2.462e+0 


1.9297e+l 


-1.020e+2 


18 


-1.2972e-2 


2.4399e-l 


-3.098e+0 


2.6888e+l 


-1.574e+2 


19 


-9.7002e-3 


1.9560e-l 


-2.643e+0 


2.4160e+l 


-1.468e+2 


20 


-2.7292e-4 


-2.2153e-2 


8.8442e-l 


-1.648e+l 


1.8126e+2 


21 


1.7454e-2 


-5.2713e-l 


1.0822e+l 


-1.530e+2 


1.4694e+3 
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TABLE 7 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR 



SNR (dB) 


I 


J 


K 


0 


-8.237e+l 


3.3306e+l 


-7.82 le+0 


1 


-7.150e+l 


3.0643e+l 


-7.946e+0 


2 


-6.261e+l 


2.8453e+l 


-8.185e+0 


3 


-5.505e+l 


2.6566e+l 


-8.540e+0 


4 


-4.838e+l 


2.4869e+l 


-9.022e+0 


5 


-3.878e+l 


2.2004e+l 


-9.490e+0 


6 


-2.857e+l 


1.8617e+l 


-1.002e+l 


7 


-1.915e+l 


1.5169e+l 


-1.070e+l 


8 


-1.053e+l 


. 1.1564e+l 


-1.150e+l 


9 


-5.058e+0 


8.997 le+0 


-1.270e+l 


10 


-2.515e+0 


7.7873e+0 


-1.450e+l 


11 


-4.7064e-l 


6.6230e+0 


-1.670e+l 


12 


2.935 le+0 


3.7027e+0 


-1.860e+l 


13 


1.0234e+l 


-4.680e+0 


- 1 ,806c+l 


14 


2.7395e+l 


-2.890e+l 


-7.690e+0 


15 


6.9656e+l 


-9.909e+l 


3.7526e+l 


16 


1.6534e+2 


-2.826e+2 


1.8421e+2 


17 


3.4745e+2 


-6.825e+2 


5.5860e+2 


18 


5.9439e+2 


-1.30e+3 


1.2141e+3 


19 


5.6268e+2 


-1.206e+3 


1.0247e+3 


20 


-1.205e+3 


4.5084e+3 


-7.352e+3 


21 


-9.182e+3 


3.3716e+4 


-5.530e+4 
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TABLE 8: EXPONENTIAL COEFFICIENTS, a = 0.948, h = 0.393, FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xjj 


A 


B 


C 


0 


.01 


5.49 


-1.8377e-3 


5.357 le-2 


-6.7308e-l 


1 


.015 


5.665 


-1.1871e-3 


3.5775e-2 


-4.6488e-l 


2 


.02 


5.865 


-7.5975e-4 


2.3744e-2 


-3.2007e-l 


3 


.025 


6.085 


-4.8463e-4 


1.5736e-2 


-2.2044e-l 


4 


.03 


6.34 


-3.0099e-4 


1.0193e-2 


-1.4896e-l 


5 


.045 


6.62 


-1.5977e-4 


5.6676e-3 


-8.6809e-2 


6 


.075 


6.935 


-7.2123e-5 


2.6945e-3 


-4.3514e-2 


7 


.13 


7.29 


-2.7062e-5 


1.0714e-3 


-1.8365e-2 


8 


.255 


7.685 


-6.7088e-6 


2.844 le-4 


-5.2386e-3 


9 


.495 


8.125 


-1.0933e-6 


4.9927e-5 


-9.9742e-4 


10 


.86 


8.625 


-1.451 le-7 


6.6944e-6 


-1.367 le-4 


11 


1.315 


9.18 


4.3950e-8 


-3.0324e-6 


8.5833e-5 


12 


1.855 


9.805 


9.5767e-8 


-6.3203e-6 


1.8027e-4 


13 


2.485 


10.5 


1.0725e-7 


-7.5820e-6 


2.3393e-4 


14 


3.21 


11.285 


1.0874e-7 


-8.368 le-6 


2.8277e-4 


15 


4.035 


12.165 


1.1220e-7 


-9.4587e-6 


3.5192e-4 


16 


4.975 


13.15 


1.1575e-7 


-1.0742e-5 


4.4164e-4 


17 


6.03 


14.255 


1.1395e-7 


-1.1682e-5 


5.3215e-4 


18 


7.22 


15.5 


1.0529e-7 


-1.1953e-5 


6.041 le-4 


19 


8.55 


16.895 


8.9424e-8 


-1.1243e-5 


6.3014e-4 


20 


10.045 


18.46 


6.6838e-8 


-9.2765e-6 


5.7415e-4 


21 


11.715 


20.215 


3.9626e-8 


-5.9821e-6 


4.0160e-4 
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TABLE 8 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393, FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


4.7703e+0 


-2.096e+l 


5.9145e+l 


-1.074e+2 


1.2262e+2 


1 


3.4094e+0 


-1.552e+l 


4.5384e+l 


-8.556e+l 


1.0166e+2 


2 


2.4362e+0 


-1.151e+l 


3.5001e+l 


-6.866e+l 


8.5004e+l 


3 


1.7443e+0 


-8.574e+0 


2.7124e+l 


-5.54 le+1 


7.1507e+l 


4 


1.2299e+0 


-6.309e+0 


2.0839e+l 


-4.446e+l 


5.9946e+l 


5 


7.5182e-l 


-4.050e+0 


1.4068e+l 


-3.163e+l 


4.5063e+l 


6 


3.9790e-l 


-2.267e+0 


8.3527e+0 


-1.999e+l 


3.0472e+l 


7 


1.7864e-l 


-1.086e+0 


4.286 le+0 


-1.105e+l 


1.8296e+l 


8 


5.5014e-2 


-3.6341e-l 


. 1.5717e+0 


-4.492e+0 


8.3925e+0 


9 


1.1476e-2 


-8.4234e-2 


4.1246e-l 


-1.366e+0 


3.0493e+0 


10 


1.6566e-3 


-1.3502e-2 


7.8663e-2 


-3.3111e-l 


9.8826e-l 


11 


-1.3158e-3 


U992e-2 


-6.6464e-2 


2.1659e-l 


-3.4775e-l 


12 


-2.9094e-3 


2.9215e-2 


-1.8910e-l 


7.9124e-l 


-2.082e+0 


13 


-4.1304e-3 


4.5996e-2 


-3.3587e-l 


1.6214e+0 


-5.087e+0 


14 


-5.5089e-3 


6.8306e-2 


-5.6153e-l 


3.0925e+0 


-1.125e+l 


15 


-7.5940e-3 


1.0504e-l 


-9.7154e-l 


6.0774e+0 


-2.538e+l 


16 


-1.0580e-2 


1.6332e-l 


-1.696e+0 


1.1991e+l 


-5.699e+l 


17 


-1.4170e-2 


2.4407e-l 


-2.840e+0 


2.2590e+l 


-1.213e+2 


18 


-1.7889e-2 


3.4355e-l 


-4.469e+0 


3.9865e+l 


-2.408e+2 


19 


-2.0724e-2 


4.4271e-l 


-6.417e+0 


6.3885e+l 


-4.314e+2 


20 


-2.0858e-2 


4.9230e-l 


-7.884e+0 


8.6732e+l 


-6.469e-f-2 


21 


-1.5768e-2 


4.0030e-l 


-6.853e-fO 


7.9908e+l 


-6.246e+2 
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TABLE 8 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393, FOR .x^ 



SNR (dB) 


I 


J 


K 


0 


-8.451e+l 


3.3738e+l 


-7.844e+0 


I 


-7.323e+l 


3.1009e+l 


-7.967e+0 


2 


-6.402e+l 


2.8766e+l 


-8.205e+0 


3 


-5.628e+l 


2.6850e+l 


-8.5595e+0 


4 


-4.939e+l 


2.5114e+l 


-9.040e+0 


5 


-3.954e+l 


2.2207e+l 


-9.508e+0 


6 


-2.903e+l 


1.8764e+l 


-1.004e+l 


7 


-1.940e+l 


1.5275e+l . 


-1.071e+l 


8 


-1.060e+l 


1.1637e+l 


-1.152e+l 


9 


-5.119e+0 


9.1053e+0 


-1.272e+l 


10 


-2.61 2e+0 


7.9441e+0 


-1.454e+l 


11 


-6.0355e-l 


6.8294e+0 


-1.675e+l 


12 


2.6179e+0 


4.1167e+0 


-1.874e+l 


13 


9.3145e+0 


-3.536e+0 


-1.858e+l 


14 


2.5322e+l 


-2.619e+l 


-9.138e+0 


15 


6.7212e+l 


-9.598e+l 


3.5977e+l 


16 


1.7376e+2 


-3.007e+2 


2.0102e+2 


17 


4.2075e+2 


-8.438e+2 


7.1562e+2 


18 


9.4237e+2 


-2.150e+3 


2.1391e+3 


19 


1.8908e+3 


-4.849e+3 


5.4842e+3 


20 


3.1296e+3 


-8.856e+3 


1.1074644 


21 


3.1156e+3 


-8.880e+3 


1.0771644 
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TABLE 9: EXPONENTIAL COEFFICIENTS, a = 0.96043, b = 0.39782, FOR 



SNR (dB) 


Range 
Lower 
Limit xl 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.015 


5.56 


-1.4228e-3 


4.2082e-2 


-5.3667e-l 


1 


.015 


5.74 


-1.0458e-3 


3.1934e-2 


-4.2043e-l 


2 


.02 


5.94 


-6.7273e-4 


2.1292e-2 


-2.9066e-l 


3 


.025 


6.165 


-4.2793e-4 


1.4076e-2 


-1.9977e-l 


4 


.035 


6.415 


-2.4809e-4 


8.5119e-3 


-1.2606e-l 


5 


.05 


6.7 


-1.3346e-4 


4.7964e-3 


-7.4453e-2 


6 


.075 


7.02 


-6.4488e-5 


2.4384e-3 


-3.9852e-2 


7 


.135 


7.38 


'-2.3286e-5 


9.3387e-4 


-1.6217e-2 


8 


.260 


7.78 


-5.869 le-6 


2.5196e-4 


-4.6999e-3 


9 


.505 


8.23 


-9.4667e-7 


4.3803e-5 


-8.8678e-4 


10 


.870 


8.73 


-1.3013e-7 


6.0822e-6 


-1.2580e-4 


11 


1.335 


9.295 


3.9229e-8 


-2.7342e-6 


7.8260e-5 


12 


1.88 


9.925 


8.4119e-8 


-5.625 le-6 


1.6255e-4 


13 


2.52 


10.635 


9.4973e-8 


-6.7975e-6 


2.1234e-4 


14 


3.255 


11.43 


9.6305e-8 


-7.5038e-6 


2.5675e-4 


15 


4.095 


12.32 


9.9680e-8 


-8.5092e-6 


3.2058e-4 


16 


5.04 


13.32 


1.0205e-7 


-9.5910e-6 


3.9933e-4 


17 


6.11 


14.44 


1.0037e-7 


-1.0422e-5 


4.8078e-4 


18 


7.315 


15.695 


9.2926e-8 


-1.0683e-5 


5.4675e-4 


19 


8.665 


17.11 


7.8618e-8 


-1.0009e-5 


5.6808e-4 


20 


10.18 


18.695 


5.8443e-8 


-8.2119e-6 


5.1454e-4 


21 


11.87 


20.475 


3.4172e-8 


-5.2174e-6 


3.5415e-4 
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TABLE 9 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.96043, b = 0.39782, FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


3.8628e+0 


-1.725e+l 


4.9533e+l 


-9.168e+l 


1.0700e+2 


I 


3.1241e+0 


-1.440e+l 


4.2684e+l 


-8.152e+l 


9.8112e+l 


2 


2.2405e+0 


-1.072e+l 


3.3008e+l 


-6.556e+l 


8.2179e+l 


3 


1.6013e+0 


-7.973e+0 


2.555 le+1 


-5.287e+l 


6.9095e+l 


4 


1.0552e+0 


-5.490e+0 


1.8407e+l 


-3.990e+l 


5.4754e+l 


5 


6.5366e-l 


-3.57 le+0 


1.2587e+l 


-2.874e+l 


4.1638e+l 


6 


3.6878e-l 


-2.127e+0 


7.9267e+0 


-1.919e+l 


2.9595e+l 


7 


1.5984e-l 


-9.8490e-l 


3.9402e+0 


-1.030e+l 


1.7310e+l 


8 


4.9987e-2 


-3.3447e-l 


1.4654e+0 


-4.244e+0 


8.0344e+0 


9 


1.0341e-2 


-7.6957e-2 


3.8215e-l 


-1.284e+0 


2.9082e+0 


10 


1.5431e-3 


-1.2717e-2 


7.4842e-2 


-3.1804e-l 


9.5835e-l 


11 


-1.2140e-3 


1.1202e-2 


-6.2904e-2 


2.0795e-l 


-3.4027e-l 


12 


-2.6576e-3 


2.7033e-2 


-1.7723e-l 


7.5112e-l 


-2.002e+0 


13 


-3.7962e-3 


4.2808e-2 


-3.1657e-l 


1.5478e+0 


-4.919e+0 


14 


-5.0650e-3 


6.3598e-2 


-5.2950e-l 


2.9536e+0 


-1.089e+l 


15 


-7.0056e-3 


9.8143e-2 


-9.1939e-l 


5.8257e+0 


-2.465e+l 


16 


-9.6877e-3 


1.5146e-l 


-1.593e+0 


U404e+1 


-5.489e+l 


17 


-1.2965e-2 


2.2616e-l 


-2.665e+0 


2.1469e+l 


-1.168e+2 


n 


-1.6396e-2 


3.1887e-l 


-4.200e+0 


3.7946e+l 


-2.32 le+2 


19 


-1.8919e-2 


4.0923e-l 


-6.006e+0 


6.0547e+l 


-4.140e+2 


20 


-1.8922e-2 


4.5206e-l 


-7.327e+0 


8.1574e+l 


-6.157e+2 


21 


-1.4054e-2 


3.6043e-l 


-6.229e+0 


7.3261e+l 


-5.768e+2 
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TABLE 9 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.96043, b = 0.39782 FOR 



SNR (dB) 


I 


J 


K 


0 


-7.579e+l 


3.1419641 


-7.680640 


I 


-7.158e+l 


3.06836-41 


-7.988640 


2 


-6.265e+l 


2.848064-1 


-8.2276-40 


3 


-5.506e+l 


2.65816-41 


-8.5826-40 


4 


-4.606e+l 


2.40786-41 


-8.9716-40 


5 


-3.725e+l 


2.14516-41 


-9.463640 


6 


-2.85 le+1 


1.8606641 


-1.0066-41 


7 


-1.865e+l 


1.4966641 


-1.0716-41 


8 


-1.029e+l 


. 1.1467641 


-1.153641 


9 


-4.960e+0 


8.9704640 


-1.2736-41 


10 


-2.555e+0 


7.8490640 


-1.4566-41 


11 


-5.788e-l 


6.7328640 


-1.6766-41 


12 


2.5493e+0 


4.0655640 


-1.8766-41 


13 


9.1301640 


-3.554640 


-1.8566-41 


14 


2.4821e+l 


-2.605641 


-9.029640 


15 


6.6140e+l 


-9.5796-41 


3.6711641 


16 


1.6950642 


-2.970642 


2.0113642 


17 


4.101364-2 


-8.3296-42 


7.1543642 


18 


9.199464-2 


-2.1266-43 


2.1417643 


19 


1.8372643 


-4.770643 


5.4618643 


20 


3.0134643 


-8.6256-43 


1.0908644 


21 


2.8923643 


-8.259643 


9.9716643 
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(x)od 



Figure 29: Third Order Curve-Fits for p^ix), with Residual Errors 
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IV. PDF FOR THE ENVELOPE DETECTION 
APPROXIMATION al/l + b\Q\ OF A TARGET TEST 

CELL 



A. CLOSED FORM PDF BY ANALYTICAL METHODS 

The envelope detection approximation given by (2) is a simplified version of 
the approximation considered so far in this thesis. It does not have the Max and Min 
operations and can therefore be calculated faster than x^. Analytically, Xg also has the 
advantage that the probability density function can be derived without the need for error 
functions. As can be recalled in the pdf development for x^, the error function terms arose 
from the expressions for the cumulative distribution functions. These were required to 
describe the statistical distribution of the Max and Min operations. 

In this situation, the derivation of the density function of Xg for a target test cell 
requires the convolution of two simple linear transformations of the random variables I/I 
and 1(21 by the multiplying coefficients a and b. Applying (44) to equations (21) and (22), 
the two terms in the convolution are 



Pa\i\^^) = 



aj^ 



^ ■ 



+ e 



\,x 

-2 (-+«»/) 



w(x) 



(83) 



and 






lx / 



+ e 



_i(f, 

2^b 



2 

■mo) ^ 



u{x) 



(84) 



These can be convolved since I and Q are assumed to be mutually independent Gaussian 
random variables. As can be recalled from Section HI, the convolution arises from the 
addition of the two independent components al/l and b\Q\. 
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The convolution integral is given by 



Pc(^) = \Pa\i\^y)Pb\Q\^^-y)^y 



or 



X 

Pc^^) = \Pa\[\^y)Pb\Q\^^-y)^y • (ss) 

0 



The integration limits reduce as shown when the unit step functions of (83) and (84) are 
removed from the integration integral. Expanding (85), the convolution integral becomes 



Pc (-^) 



1 

2abn 







^ r2}?/n^6 + IxrriQa 


1 




(— 

o* 

L_ 


0 







~2ymjb-^2ymQcr 


1 




ab 


b _ 


+ . J 






J 



X 



2a^b^ 



(-2xya^ 



hyW- 



2 2.2 . -2 2 
■m,a b +JC a - 






2 m 



b- 



m^Qab^ A 



IxniQab' 



lyrnfib) 



dy 



or 



Pc 



1 

2abn 



X 



1 






[— (a^ + fe^))’^+ {Amfib^ •*-2x(^ — 2mQab — 2miab^)y-¥ (^xniQab — 



mja^b^ — — m^a^b^ — 2xmQa^b ) ] 



+ 



^ ^ [- (a^ + b^) y^ + {Am,ab^ + 4n>Qa^b + 2xa^ -2mQC^b — 2miab^) y+ (-m^a^b^- 

^e2ab 

0 



X (^-tr?Q(^b^-2xmQC^b ) ] 



dy + 
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X — [- + {2xa - Inioab - Irrijab^) y + {AxmQab-m]ab ^ - 

0 



X c^-m^QC^b^—lxmQC^b ) ] 



rf}’ + 



i ^ [— (a^ + b^)y^+ (AmQab^ + 2xd^-2mQa^b — 2miab^)y+ (—rnifa^b^- 
0 

X c?-rr?Q<?l?—2xmQC?b ) ] 



Each of the four integral terms in (86) are now of the form 



( 86 ) 






+ By + C 



dy 



( 87 ) 



which integrate in closed form to 

^ [ACA-^ 

Ju 4^ J 

2j^ 

Applied to (86), this gives the result 



f ff ^ \ ^ \ \ 



( 88 ) 



Pt(^) = 



1 


Jn 

p 


~AC^A-BA 


L 4A J 


2abn 


2^A 





‘ ^2^] A 2j=A 



} + 



Jn 

Jn 

2]^' 



AC^-{By 

AA 



AC,A-Bp 



L AA 



f \ (2Ax + B^\ 

27^ j’ 
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Jn 



'AC^A-B^ 
_ 44 



{erf 



f Ba \ 



(2Ax + 

-erf\—T=J^\) 



\2j^J { 2j^ 



(89) 



where 



A = -ia^ + b^) 

2a r 



(90) 



5 j = -\-2xa^ -2mQO'b -2mjal?') { — ^) 

2a b 



(91) 



^2 = {Arrijal?' AniQa'b ■^ 2 x 0 ' -2mQO'b -2mjal?') { — ^) (92) 

2a b 



= (2xa^ -2mQa^b -2m.ab^) { — ^) 

2a b^ 



( 93 ) 



£4 = {AniQab^ + 2x0' -2mQO'b-2mjab^) ( — ^) 

2a b 



(94) 



^ 2 . 2 2.2 ..2 2 2 2.2 ... 2 . w 1 \ 

Cj = {AxniQa b — nija b -x a - niQa b -2xniQa b ) ( — ^) ^ 95 ^ 

2a b 



^ , 2 2.2 .2 2 2 2.2 2. ^ 1 ^ 

C 2 = {-m^a b -X a -rriQa b -2xmQa b ) ( — ^) 

2a b 



(96) 



C 3 = {AxmQO'b—ri^(P'l?'—^(^—iri^a'lP'—2xmQa'b ){ — ^— ^) (97) 

2a b 



^ , 2 2.2 .2 2 2 2,2 O ' 2 . . . 1 . 

C 4 = {-nija b -X a -niQa b -2xmQa b ) ( — ^) 

2a^b^ 



(98) 
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B. CLOSED FORM PDF BY CURVE-FITTING 

Using the same procedure as that used for Section IH.B, the density function (x) 
is also obtained by numerical convolution for a number of signal-to-noise ratios. The results 
are curve-fitted using an exponential expansion of the form of (81). The exact coefficient 
values for each SNR are in Tables 10 to 16 in order to provide a quick and easy reference. 
This procedure allows the two envelope approximations and to be compared in terms 
of their representing coefficients. 
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TABLE 10: EXPONENTIAL COEFFICIENTS, a = 1.0, Z? = 1.0 FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.005 


7.48 


-5.6415e-5 


2.2342e-3 


-3.8108e-2 


1 


.005 


7.7 


-4.2441e-5 


1.7303e-3 


-3.0380e-2 


2 


.005 


7.95 


-3.1021e-5 


1.3057e-3 


-2.3666e-2 


3 


.005 


8.23 


-2.2070e-5 


9.6143e-4 


-1.8035e-2 


4 


.005 


8.55 


-1.5122e-5 


6.8408e-4 


-1.3325e-2 


5 


.005 


8.915 


-9.9361e-6 


4.6842e-4 


-9.5072e-3 


6 


.01 


9.325 


-5.5082e-6 


2.7178e-4 


-5.7740e-3 


7 


.025 


9.785 


-2.532 le-6 


1.3135e-4 


-2.935 le-3 


8 


.09 


10.305 


-6.6428e-7 


3.6566e-5 


-8.6903e-4 


9 


.345 


10.895 


-2.3983e-8 


1.488 le-6 


-4.1031e-5 


10 


.835 


11.560 


2.2540e-8 


-1.3770e-6 


3.5360e-5 


11 


1.425 


12.305 


2.7403e-8 


-1.8370e-6 


5.2550e-5 


12 


2.09 


13.145 


2.7140e-8 


-1.9782e-6 


6.1884e-5 


13 


2.835 


14.090 


2.3047e-8 


-1.8092e-6 


6.0928e-5 


14 


3.655 


15.155 


1.6102e-8 


-1.3247e-6 


4.6002e-5 


15 


4.57 


16.355 


7.2909e-9 


-5.3355e-7 


1.3083e-5 


16 


5.585 


17.700 


-1.2846e-9 


4.0317e-7 


-3.4795e-5 


17 


6.715 


19.210 


-7.9937e-9 


1.2949e-6 


-9.0238e-5 


18 


7.97 


20.910 


-1.1564e-8 


1.9114e-6 


-1.3899e-4 


19 


9.375 


22.820 


-1.2178e-8 


2.1652e-6 


-1.7052e-4 


20 


10.945 


24.970 


-1.0577e-8 


2.0534e-6 


-1.7722e-4 


21 


12.700 


27.380 


-7.9595e-9 


1.6944e-6 


-1.6077e-4 
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TABLE 10 (cont’d): EXPONENTIAL COEFFICIENTS, a=l.0,b = 1.0 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


3.6625e-l 


-2.179e+0 


8.3030e+0 


-2.0276+1 


3.08856+1 


1 


3.0056e-l 


-1.841e+0 


7.2174e+0 


-1.8136+1 


2.83776+1 


2 


2.4169e-l 


-1.528e+0 


6.1808640 


-1.6016+1 


2.58086+1 


3 


1.9059e-l 


-1.246e+0 


5.21486+0 


-1.3966+1 


2.32206+1 


4 


1.4619e-l 


-9.9219e-l 


4.30706+0 


-1.1956+1 


2.05736+1 


5 


1.0867e-l 


-7.6809e-l 


3.47076+0 


-l.OOle+1 


1.79046+1 


6 


6.9094e-2 


-5.1141e-l 


2.42046+0 


-7.3156+0 


1.37036+1 


7 


3.6968e-2 


-2.8825e-l 


1.43856+0 


-4.5906+0 


9.0966e+0 


8 


1.1674e-2 


-9.7447e-2 


■ 5.22936-1 


-1.8046+0 


3.89726+0 


9 


6.5957e-4 


-6.7715e-3 


4.54726-2 


-1.97116-1 


5.45686-1 


10 


-4.9293e-4 


4.0194e-3 


-1.93136-2 


5.27586-2 


-6.1962e-2 


11 


-8.362 le-4 


8.1102e-3 


-4.9748e-2 


1.96316-1 


-4.90756-1 


12 


-1.0857e-3 


1.1741e-2 


-8.14106-2 


3.67786-1 


-1.0726+0 


13 


-1.1490e-3 


1.3298e-2 


-9.76656-2 


4.56596-1 


-1.3186+0 


14 


-8.6775e-4 


9.4267e-3 


-5.5206e-2 


9.8004e-2 


7.9852e-l 


15 


-1.6755e-5 


-5.7953e-3 


1.39296-1 


-1.6576+0 


1.15716+1 


16 


1.4906e-3 


-3.7844e-2 


6.1262e-l 


-6.494e+0 


4.48986+1 


17 


3.5860e-3 


-9.0311e-2 


1.50926+0 


-1.697e+l 


1.2689e+2 


n 


5.8579e3 


-1.5848e-l 


2.87556+0 


-3.5436+1 


2.92716+2 


19 


7.8313e-3 


-2.3215e-l 


4.63946+0 


-6.3276+1 


5.81386+2 


20 


8.9499e-3 


-2.9274e-l 


6.47646+0 


-9.8106+1 


1.00446+3 


21 


8.9490e-3 


-3.2342e-l 


7.92586+0 


-1.3336+2 


1.51956+3 
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TABLE 10 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, = 1.0 FOR 



SNR (dB) 


I 


J 


K 


0 


-2.809e+l 


1.468 le+1 


-5.179e-h0 


1 


-2.648e+l 


1.4288e+l 


-5.414e+0 


2 


-2.476e+l 


1.3863e+l 


-5.713e+0 


3 


-2.295e+l 


1.3408e+l 


-6.093e+0 


4 


-2.102e+l 


1.2915e+l 


-6.576e+0 


5 


-1.898e+l 


1.2390e+l 


-7.189e+0 


6 


-1.531e+l 


1.1092e+l 


-7.850e+0 


7 


-1.092e+l 


9.3509e+0 


-8.633e+0 


8 


-5.396e+0 


6.8162e+0 


-9.493e+0 


9 


-1.347e+0 


4.7240e+0 


-1.065e+l 


10 


-4.6656e-l 


4.4785e+0 


-1.259e+l 


11 


3.1129e-l 


4.1940C40 


-1.494e+l 


12 


1.4999e+0 


3.4208e+0 


-1.768e+l 


13 


1.7129e+0 


4.4670e+0 


-2.264e+l 


14 


-6.206e+0 


2.1808e+l 


-4.270e+l 


15 


-4.852e+l 


1.1780e+2 


-1.421e+2 


16 


-1.954e+2 


4.908 le+2 


-5.612e+2 


17 


-6.054e+2 


1.6719e+3 


-2.056e+3 


18 


-1.553e+3 


4.7866e+3 


-6.542e+3 


19 


-3.445e+3 


1.1896e44 


-1.822e+4 


20 


-6.65 le+3 


2.5735e+4 


-4.423e+4 


21 


-1.123e+4 


4.8574e+4 


-9.350e+4 
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TABLE 11: EXPONENTIAL COEFFICIENTS, a = l.O,b = 0.5 FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.005 


6.11 


-4.0739e-4 


1.3173e-2 


-1.8349e-l 


1 


.005 


6.295 


-3.0346e-4 


1.0113e-2 


-1.4517e-l 


2 


.005 


6.5 


-2.2156e-4 


7.6256e-3 


-1.1306e-l 


3 


.005 


6.73 


-1.5757e-4 


5.6153e-3 


-8.6196e-2 


4 


.005 


6.99 


-1.0845e-4 


4.0133e-3 


-6.3966e-2 


5 


.005 


7.285 


-7.1773e-5 


2.7669e-3 


-4.5935e-2 


6 


.005 


7.615 


-4.5710e-5 


1.8409e-3 


-3.1923e-2 


7 


.015 


7.99 


-2.1107e-5 


8.9384e-4 


-1.6306e-2 


8 


.045 


8.41 


-7.1578e-6 


3.2085e-4 


-6.2039e-3 


9 


.2 


8.88 


-5.5374e-7 


2.7465e-5 


-5.9376e-4 


10 


.605 


9.41 


-1.1347e-7 


6.4902e-6 


-1.6294e-4 


11 


1.14 


10.01 


-9.2300e-8 


5.475 le-6 


-1.4238e-4 


12 


1.775 


10.68 


-1.1528e-8 


7.3671e-7 


-2.0688e-5 


13 


2.505 


11.435 


6.037e-ll 


-1.5928e-8 


8.5076e-7 


14 


3.32 


12.285 


4.216e-10 


-4.5317e-8 


1.9599e-6 


15 


4.235 


13.235 


1.6105e-9 


-1.5368e-7 


6.3853e-6 


16 


5.26 


14.305 


3.3336e-9 


-3.3600e-7 


1.4918e-5 


17 


6.405 


15.51 


5.4627e-9 


-5.9855e-7 


2.9026e-5 


18 


7.685 


16.86 


7.5495e-9 


-9.0561e-7 


4.8190e-5 


19 


9.11 


18.38 


8.5007e-9 


-1.1104e-6 


6.4304e-5 


20 


10.7 


20.085 


6.5009e-9 


-8.8275e-7 


5.2265e-5 


21 


12.47 


22.0 


-8.253e-10 


3.4981e-7 


-4.3519e-5 
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TABLE 11 (cont’d): EXPONENTIAL COEFFICIENTS, a = l.0,b = 0.5 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


1.4406e+0 


-7.0065e+0 


2.1857e+l 


-4.382e+l 


5.5080e+l 


1 


1.1748e-fO 


-5.889e+0 


1.8933e+l 


-3.910e+l 


5.0587e+l 


2 


9.4483e-l 


-4.891e+0 


1.6235e+l 


-3.460e+l 


4.6122e+l 


3 


7.4575e-l 


-3.996e+0 


1.3724e+l 


-3.024e+l 


4.1612e+l 


4 


5.7453e-l 


-3.195e+0 


1.1383e+l 


-2.60e+l 


3.7020e+l 


5 


4.2966e-l 


•2.488c+0 


9.2214e+0 


-2.189e+l 


3.2366e+l 


6 


3.1182e-l 


-1.885e+0 


7.2899e+0 


-1.805e+l 


2.7784e+l 


7 


1.6768e-l 


-1.068e+0 


4.3592e+0 


-1.141e+l 


1.8652e+l 


8 


6.7743e-2 


-4.5935e-l 


2.0043e+0 


-5.643e+0 


1.0016e+l 


9 


7.3490e-3 


-5.7525e-2 


2.9685e-l 


-1.020e+O 


2.3035e+0 


10 


2.3593e-3 


-2.1763e-2 


1.3319e-l 


-5.4507e-l 


1.4657e+0 


11 


2.1334e-3 


-2.0347e-2 


1.2871e-l 


-5.4448e-l 


1.5139e+0 


12 


3.3607e-4 


-3.4953e-3 


2.4246e-2 


-1.1255e-l 


3.4086e-l 


13 


-2.0985e-5 


2.9302e-4 


-2.562 le-3 


1.5345e-2 


-6.6124e-2 


14 


-4.5959e-5 


6.5555e-4 


-6.0484e-3 


3.7806e-2 


-1.6376e-l 


15 


-1.5180e-4 


2.2859e-3 


-2.2859e-2 


1.5515e-l 


-7.1334e-l 


16 


-3.8360e-4 


6.3234e-3 


-6.9900e-2 


5.2649e-l 


-2.679e+0 


17 


-8.1963e-4 


1.4920e-2 


-1.8304e-l 


1.5350e+0 


-8.705e+0 


18 


-1.4970e-3 


3.0057e-2 


-4.0764e-l 


3.7849e+0 


-2.378e+l 


19 


-2.1720e-3 


4.7343e-2 


-6.9534e-l 


6.9659e+0 


-4.698e+l 


20 


-1.7601e-3 


3.6752e-2 


-4.8213e-l 


3.7350e+0 


-1.254e+l 


21 


2.7514e-3 


-1.0509e-l 


2.6089e+0 


-4.326e+l 


4.7734e+2 
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TABLE 11 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5 FOR 



SNR (dB) 


I 


J 


K 


0 


-4.150e+l 


1.7671e+l 


-4.652e+0 


1 


-3.914e+l 


1.7202e+l 


-4.887e+0 


2 


-3.669e+l 


1.6707e+l 


-5.187e+0 


3 


-3.410e+l 


1.6176e+l 


-5.568e+0 


4 


-3.134e+l 


1.5602e+l 


-6.053e+0 


5 


-2.842e+l 


1.4990e+l 


-6.668e+0 


6 


-2.541e+l 


1.4362e+l 


-7.449e+0 


7 


-1.835e+l 


1.2085e+l 


-8.220e+0 


8 


-1.106e+l 


9.4920e+0 


-9.173e+0 


9 


-3.706e+0 


6.4565e+0 


-1.029e+l 


10 


-2.902e+0 


6.5895e+0 


-1.236e+l 


11 


-3.069e+0 


7.3277e+0 


-1.507e+l 


12 


-1.069e+0 


5.9483e+0 


-1.750e+l 


13 


-2.3636e-l 


5.5735e+0 


-2.103e+l 


14 


4.4512e-2 


5.7907e+0 


-2.576e+l 


15 


1.7067e+0 


3.6425e+0 


-2.988e+l 


16 


8.401 le+0 


-8.736e+0 


-2.633e+l 


17 


3.1562e+l 


-5.952e+l 


1.4733e+l 


18 


9.6423e+l 


-2.206e+2 


1.8119e+2 


19 


2.035 le+2 


-5.022e+2 


4.8969e+2 


20 


-3.550e+l 


4.5997e+2 


-1.221e+3 


21 


-3.373e+3 


1.3841e+4 


-2.518e+4 
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TABLE 12: EXPONENTIAL COEFFICIENTS, a=l.0,b = 0.25 FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit Xy 


A 


B 


C 


0 


.005 


5.665 


-8.8816e-4 


2.6767e-2 


-3.4782e-l 


1 


.005 


5.835 


-6.6626e-4 


2.0687e-2 


-2.7696e-l 


2 


.005 


6.025 


-4.8777e-4 


1.5641e-2 


-2.1623e-l 


3 


.005 


6.24 


-3.4597e-4 


1.1489e-2 


-1.6450e-l 


4 


.005 


6.485 


-2.3601e-4 


8.1444e-3 


-1.2116e-l 


5 


.005 


6.755 


-1.5595e-4 


5.6048e-3 


-8.6829e-2 


6 


.005 


7.06 


-9.8414e-5 


3.6966e-3 


-5.985 le-2 


7 


.01 


7.405 


-5.065 le-5 


2.0015e-3 


-3.4111e-2 


8 


.025 


7.79 


-2.1057e-5 


8.8247e-4 


-1.5975e-2 


9 


.105 


8.23 


-5.4285e-6 


2.4762e-4 


-4.9015e-3 


10 


.375 


8.715 


-3.3928e-6 


1.6444e-4 


-3.4603e-3 


11 


.795 


9.27 


-4.7802e-7 


2.4862e-5 


-5.6183e-4 


12 


1.38 


9.885 


4.4594e-8 


-2.639 le-6 


6.8575e-5 


13 


2.04 


10.58 


-1.521e-10 


5.1163e-9 


4.0578e-8 


14 


2.775 


11.36 


-1.0876e-9 


7.1638e-8 


-1.9784e-6 


15 


3.6 


12.24 


-2.5077e-9 


1.9182e-7 


-6.3354e-6 


16 


4.525 


13.225 


-5.6096e-9 


4.9085e-7 


-1.8818e-5 


17 


5.55 


14.335 


-1.2042e-8 


1.1958e-6 


-5.2465e-5 


18 


6.7 


15.580 


1 -2.441e-8 


2.7380e-6 


-1.3639e-4 


19 


7.985 


16.975 


-4.6017e-8 


5.8150e-6 


-3.2740e-4 


20 


9.425 


18.545 


-7.9046e-8 


1.1236e-5 


-7.1324e-4 


21 


11.03 


20.310 


-1.2136e-7 


1.9393e-5 


-1.3860e-3 
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TABLE 12 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, Z? = 0.25 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


2.5507e+0 


-1.160e+l 


3.3900e+l 


-6.368e+l 


7.4832e+l 


1 


2.0924e+0 


-9.806e+0 


2.9504e+l 


-5.705e-Hl 


6.8917e+l 


2 


1.6869e+0 


-8.161e+0 


2.5344e+l 


-5.054e+l 


6.289 le+1 


3 


1.3288e+0 


-6.655e+0 


2.1388e+l 


-4.41 le+1 


5.6686e+l 


4 


1.0168e+0 


-5.289e+0 


1.7647e+l 


-3.776e+l 


5.0307e+l 


5 


7.5875e-l 


-4.109e+0 


1.4270e+l 


-3.177e+l 


4.4015e+l 


6 


5.4657e-l 


-3.094e+0 


1.1226e+l 


-2.612e+l 


3.7817e+l 


7 


3.2817e-l 


-1.959e+0 


7.5072e+0 


-1.849e+l 


2.8422e+l 


8 


1.6358e-l 


-1.042e+0 


4.2776e+0 


-1.134e+l 


1.8905e+l 


9 


5.5182e-2 


-3.8912e-l 


1.7836e+0 


-5.336e+0 


1.0174e+l 


10 


4.1443e-2 


-3.1117e-l 


1.5205e+0 


-4.858e+0 


9.9220e+0 


11 


7.2312e-3 


-5.8378e-2 


3.0678e-l 


-1.054e+0 


2.3062e+0 


12 


-1.0281e-3 


9.8245e-3 


-6.2324e-2 


2.647 le-1 


-7.4149e-l 


13 


-4.7141e-6 


1.0104e-4 


-1.0702e-3 


6.0659e-3 


-1.8244e-2 


14 


2.9113e-5 


-2.3517e-4 


8.9522e-4 


-1.5983e-4 


-9.0667e-3 


15 


1.1757e-4 


-1.3360e-3 


9.5165e-3 


-4.2298e-2 


1.1636e-l 


16 


4.1422e-4 


-5.765 le-3 


5.2732e-2 


-3.2019e-l 


1.2791e+0 


17 


1.3361e-3 


-2.1818e-2 


2.3822e-l 


-1.759e+0 


8.6865e+0 


18 


3.9680e-3 


-7.4579e-2 


9.4519e-l 


-8.175e+0 


4.7653e+l 


19 


1.0808e-2 


-2.3151e-l 


3.3602e+0 


-3.345e+l 


2.2555e+2 


20 


2.6613e-2 


-6.4615e-l 


1.0663e+l 


-1.211e+2 


9.3400e+2 


21 


5.8329e-2 


-1.600e+0 


2.9900e+l 


-3.852e+2 


3.3789e+3 
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TABLE 12 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.25 FOR 



SNR (dB) 


I 


J 


K 


0 


-5.1980e+l 


1.9376641 


-4.037640 


1 


-4.9116e+l 


1.8898641 


-4.274640 


2 


-4.6088e+l 


1.8394641 


-4.575640 


3 


-4.2853e+l 


1.7858641 


-4.960640 


4 


-3.9406e+l 


1.7296641 


-5.448640 


5 


-3.5882e+l 


1.6743641 


-6.069640 


6 


-3.2295e+l 


1.6216641 


-6.859640 


7 


-2.6005e+l 


1.4709641 


-7.734640 


8 


-1.9085e+l 


1.2864e4l 


-8.804640 


9 


-1.2114e+l 


1.0796641 


-1.012641 


10 


-1.2688e+l 


1.2192641 


-1.238641 


11 


-3.4700e+0 


6.6301640 


-1.305641 


12 


8.6893e-l 


3.7169640 


-1.473641 


13 


-4.0848e-l 


5.5502640 


-1.873641 


14 


-4.1072e-l 


6.1537640 


-2.297641 


15 


-6.2911e-l 


7.0375640 


-2.830641 


16 


-3.6775e+0 


1.2364641 


-3.780641 


17 


-2.7938e+l 


5.9044641 


-8.393641 


18 


- 1.797 le+2 


4.0290642 


-4.359642 


19 


-9.8619e+2 


2.5331643 


-2.937643 


20 


-4.6847e+3 


1.3805644 


-1.820644 


21 


-1.9309644 


6.4921644 


-9.759644 
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TABLE 13: EXPONENTIAL COEFFICIENTS, a = 1.0, 6 = 3/8 FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.005 


5.865 


-5.9627e-4 


1.8544e-2 


-2.4855e-l 


1 


.005 


6.04 


-4.467 le-4 


1.4312e-2 


-1.9764e-l 


2 


.005 


6.24 


-3.248 le-4 


1.0754e-2 


-1.5346e-l 


3 


.005 


6.46 


-2.3 13 le-4 


7.9286e-3 


-1.1713e-l 


4 


.005 


6.71 


-1.5888e-4 


5.6558e-3 


-8.6757e-2 


5 


.005 


6.99 


-1.0527e-4 


3.9019e-3 


-6.2319e-2 


6 


.005 


7.31 


-6.6350e-5 


2.5708e-3 


-4.2912e-2 


7 


.01 


7.665 


-3.4146e-5 


1.3895e-3 


-2.4367e-2 


8 


.035 


8.065 


-1.1382e-5 


4.9132e-4 


-9.1576e-3 


9 


.155 


8.515 


-1.1236e-6 


5.4691e-5 


-1.1618e-3 


10 


.505 


9.02 


-6.5642e-7 


3.4302e-5 


-7.8293e-4 


11 


1 


9.59 


-3.0786e-7 


1.7138e-5 


-4.1691e-4 


12 


1.61 


10.23 


4.6154e-9 


-3.0650e-7 


8.9859e-6 


13 


2.315 


10.95 


1.1272e-9 


-7.795 le-8 


2.3756e-6 


14 


3.105 


11.76 


-2.835e-ll 


1.9028e-9 


-6.2129e-8 


15 


3.99 


12.67 


-4.250e-ll 


2.9915e-9 


-9.9091e-8 


16 


4.985 


13.69 


-3.001e-ll 


1.7139e-9 


-3.9367e-8 


17 


6.095 


14.84 


3.854e-ll 


-5.8870e-9 


3.4002e-7 


18 


7.345 


16.125 


2.262e-10 


-2.938 le-8 


1.6541e-6 


19 


8.74 


17.575 


6.477e-10 


-8.8982e-8 


5.401 le-6 


20 


10.31 


19.2 


1.4930e-9 


-2.2408e-7 


1.4962e-5 


21 


12.06 


21.03 


3.0399e-9 


-5.0390e-7 


3.7284e-5 
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TABLE 13 (cont’d): EXPONENTIAL COEFFICIENTS, a = I.O, /> = 3/8 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


1.8793e+0 


-8.81116+0 


2.65346+1 


-5.1436+1 


6.26036+1 


1 


1. 5396640 


-7.4376+0 


2.30706+1 


-4.6046+1 


5.76306+1 


2 


1.2352640 


-6.1656+0 


1.97526+1 


-4.0696+1 


5.24946+1 


3 


9.75896-1 


-5.0406+0 


1.67066+1 


-3.5576+1 


4.73626+1 


4 


7.50476-1 


-4.0236+0 


1.38336+1 


-3.0536+1 


4.20786+1 


5 


5.61166-1 


-3.1316+0 


1.11976+1 


-2.5696+1 


3.67586+1 


6 


4.03786-1 


-2.3536+0 


8.78986+0 


-2.1056+1 


3.14206+1 


7 


2.41016-1 


-1.4786+0 


5.81146+0 


-1.4686+1 


2.31796+1 


8 


9.65166-2 


-6.32776-1 


2.67536+0 


-7.3196+0 


1.26616+1 


9 


1.41396-2 


-1.08776-1 


5.50386-1 


-1.8456+0 


4.02206+0 


10 


1.02466-2 


-8.48196-2 


4.62156-1 


-1.6706+0 


3.92486+0 


11 


5.81856-3 


-5.14146-2 


2.99376-1 


-1.1586+0 


2.91806+0 


12 


-1.52796-4 


1.66376-3 


-1.20986-2 


5.96866-2 


-1.99866-1 


13 


-4.17646-5 


4.65136-4 


-3.40656-3 


1.68076-2 


-5.86366-2 


14 


1.54396-6 


-3.21276-5 


4.74626-4 


-4.21496-3 


1.90216-2 


15 


2.36656-6 


-4.78546-5 


7.36116-4 


-7.33156-3 


4.15606-2 


16 


8.47466-7 


-2.81986-5 


6.84536-4 


-9.33006-3 


6.95596-2 


17 


-1.01176-5 


1.69396-4 


-1.55076-3 


5.77326-2 


1.79416-2 


18 


-5.29626-5 


1.06266-3 


-1.38726-2 


1.18526-1 


-6.51256-1 


19 


-1.9039e-4 


4.30826-3 


-6.52846-2 


6.70286-1 


-4.6056+0 


20 


-5.84776-4 


1.48056-2 


-2.53586-1 


2.97506+0 


-2.3616+1 


21 


-1.62106-3 


4.58486-2 


-8.81296-1 


1.16596+1 


-1.0486+2 
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TABLE 13 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8 FOR 



SNR (dB) 


I 


J 


K 


0 


-4.558e+l 


1.8432e+l 


-4.401e+0 


I 


-4.305e+l 


1.7960e+l 


-4.637e+0 


2 


-4.033e+l 


1.7446e+l 


-4.938e+0 


3 


-3.749e+l 


1.6903e+l 


-5.320e+0 


4 


-3.443e+l 


1.6317e+l 


-5.806e+0 


5 


-3.123e+l 


1.5703e+l 


-6.422e+0 


6 


-2.788e+l 


1.5068e+l 


-7.205e+0 


7 


-2.197e+l 


1.3440e+l 


-8.069e+0 


8 


-1.355e+l 


> 1.0656e+l 


-9.01 6e+0 


9 


-5.857e+0 


7.7978e+0 


-1.020e+l 


10 


-6.144e+0 


8.7885e+0 


-1.246e+l 


11 


-4.996e+0 


8.6299e+0 


-1.482e+l 


12 


1.0981e-4 


4.6520e+0 


-1.614e+l 


13 


-2.8703e-l 


5.5637e+0 


-2.01 le+1 


14 


-4.6302e-l 


6.4305e+0 


-2.497e+l 


15 


-5.4542e-l 


7.2962e+0 


-3.097e+l 


16 


-6.9732e-l 


8.4752e+0 


-3.868e+l 


17 


-6.6461e-l 


9.6411e+0 


-4.835e+l 


18 


1.7867e+0 


5.6375e+0 


-5.539e+l 


19 


1.9884e+l 


-4.11 le+1 


-1.383e+l 


20 


1.2096e+2 


-3.533e+2 


3.9927e+2 


21 


6.1282e+2 


-2.095e+3 


3.1247e+3 
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TABLE 14: EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8 FOR 



SNR (dB) 


Range 
Lower 
Limit xl 


Range 
Upper 
Limit xjj 


A 


B 


C 


0 


.005 


5.715 


-7.7098e-4 


2.3358e-2 


-3.0497e-l 


I 


.005 


5.88 


-5.8291e-4 


1.8176e-2 


-2.4427e-l 


2 


.005 


6.075 


-4.2363e-4 


1.3651e-2 


-1.8959e-l 


3 


.005 


6.29 


-3.0135e-4 


1.0055e-2 


-1.4458e-l 


4 


.005 


6.53 


-2.0812e-4 


7.2077e-3 


-1.0756e-l 


5 


.005 


6.805 


-1.3746e-4 


4.959 le-3 


-7.7080e-2 


6 


.005 


7.115 


-8.6894e-5 


3.2759e-3 


-5.3204e-2 


7 


.01 


7.46 


-4.4683e-5 


1.7690e-3 


-3.0179e-2 


8 


.035 


7.85 


-1.4787e-5 


6.2088e-4 


-1.1256e-2 


9 


.15 


8.29 


-1.4578e-6 


6.8776e-5 


-1.4159e-3 


10 


.495 


8.78 


-7.4567e-7 


3.8116e-5 


-8.5150e-4 


II 


.975 


9.335 


-3.8977e-7 


2.1162e-5 


-5.022 le-4 


12 


1.57 


9.96 


1.359e-10 


-3.4947e-8 


1.7391e-6 


13 


2.255 


10.66 


1.5245e-9 


-1.0294e-7 


3.0647e-6 


14 


3.025 


11.445 


-1.539e-10 


1.0772e-8 


-3.4041e-7 


15 


3.885 


12.330 


-1.856e-10 


1.4068e-8 


-4.8170e-7 


16 


4.855 


13.325 


-1.652e-10 


1.3146e-8 


-4.7352e-7 


17 


5.935 


14.44 


-3.901e-ll 


9.333e-10 


6.6338e-8 


18 


7.15 


15.695 


3.066e-10 


-3.9335e-8 


2.1770e-6 


19 


8.515 


17.1 


1.0559e-9 


-1.4034e-7 


8.2474e-6 


20 


10.04 


18.685 


2.502 le-9 


-3.6314e-7 


2.3456e-5 


21 


11.745 


20.46 


5.0595e-9 


-8.1189e-7 


5.8164e-5 
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TABLE 14 (conl’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


2.2460e+0 


-1.026e+l 


3.0078e+l 


-5.678e+l 


6.7297e+l 


1 


1.8517e+0 


-8.703e+0 


2.6268e+l 


-5.100e+l 


6.2110e+l 


2 


1.485 le+0 


-7.213e+0 


2.2486e+l 


-4.507e+l 


5.6570e+l 


3 


1.1725e+0 


-5.894e+0 


1.9009e+l 


-3.939e+l 


5.1028e+l 


4 


9.0511e-l 


-4.720e+0 


1.5783e+l 


-3.388e+l 


4.5409e+l 


5 


6.7543e-l 


-3.667e+0 


1.2759e+l 


-2.848e+l 


3.9637e+l 


6 


4.8705e-l 


-2.761e+0 


1.0032e+l 


-2.336e+l 


3.3911e+l 


7 


2.9035e-l 


-1.73 le+0 


6.6218e+0 


-1.626e+l 


2.4961e+l 


8 


1.1537e-l 


-7.3540e-l 


3.0225e+0 


-8.037e+0 


1.3510e+l 


9 


1.6699e-2 


-1.2450e-l 


6.1069e-l 


-1.986e+0 


4.2050e+0 


10 


1.0914e-2 


-8.8563e-2 


4.7347e-l 


-1.680e+0 


3.8842e+0 


11 


6.8405e-3 


-5.9020e-2 


3.3577e-l 


-1.270e+0 


3.1320e+0 


12 


-4.0776e-5 


5.5432e-4 


-4.7573e-3 


2.6969e-2 


-1.0335e-l 


13 


-5.2649e-5 


5.7337e-4 


-4.1166e-3 


2.0040e-2 


-6.9377e-2 


14 


6.7025e-6 


-9.4947e-5 


9.8840e-4 


-6.8804e-3 


2.6616e-2 


15 


1.0190e-5 


-1.5284e-4 


1.6856e-3 


-1.2917e-2 


6.0961e-2 


16 


1.0668e-5 


-1.7469e-4 


2.1650e-3 


-1.9178e-2 


1.0926e-l 


17 


-3.4005e-6 


5.7572e-5 


-2.7059e-4 


-3.7666e-3 


5.8466e-2 


18 


-6.837 le-5 


1.3466e-3 


-1.7348e-2 


1.4813e-l 


-8.3369e-l 


19 


-2.8181e-4 


6.1941e-3 


-9.1463e-2 


9.1942e-l 


-6.226e+0 


20 


-8.8748e-4 


2.1772e-2 


-3.6181e-l 


4.1260e+0 


-3.191e+l 


21 


-2.4490e-3 


6.7105e-2 


-1.250e+0 


1.6043e+l 


-1.400e+2 
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TABLE 14 (cont’d): EXPONENTIAL COEFFICffiNTS, a = 31/32, b = 3/8 FOR 



SNR(dB) 


I 


J 


K 


0 


-4.774e+l 


1.8856e+l 


-4.389e+0 


1 


-4.517e+l 


1.8389e+l 


-4.626e+0 


2. 


-4.231e+l 


1.7861e+l 


-4.927e+0 


3 


-3.932e+l 


1.7303e+l 


-5.309e+0 


4 


-3.615e+l 


1.6710e+l 


-5.795e+0 


5 


-3.277e+l 


1.6074e+l 


-6.4111e+0 


6 


-2.927e+l 


1.5423e+l 


-7.194e+0 


7 


-2.300e+l 


1.3733e+l 


-8.056e+0 


8 


-1.407e+l 


1.0839e+l 


-8.998e+0 


9 


-5.980e+0 


7.9006e+0 


-1.018e+l 


10 


-6.029e+0 


8.7579e+0 


-1.241e+l 


11 


-5.260e+0 


8.9190e+0 


-1.488e+l 


12 


-2.0839e-l 


5.0254e+0 


-1.629e+l 


13 


-2.9033e-l 


5.7330e+0 


-2.018e+l 


.14 


-4.9947e-l 


6.6552e+0 


-2.509e+l 


15 


-6.0647e-l 


7.5690e+0 


-3.115e+l 


16 


-8.0762e-l 


8.8203e+0 


-3.891e+l 


17 


-7.5844e-l 


9.8783e+0 


-4.840e+l 


18 


2.5662e+0 


3.7327e+0 


-5.296e+l 


19 


2.6838e+l 


-5.884e+l 


7.1317e+0 


20 


1.5986e+2 


-4.606e+2 


5.3266e+2 


21 


7.958 le+2 


-2.650e+3 


3.8734e+3 
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TABLE 15: EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393 FOR 



SNR (dB) 


Range 
Lower 
Limit xl 


Range 
Upper 
Limit Xy 


A 


B 


C 


0 


.005 


5.645 


-8.744e-4 


2.6152e-2 


-3.370e-l 


1 


.005 


5.81 


-6.586e-4 


2.028 le-2 


-2.691 le-1 


2 


.005 


6.0 


-4.807e-4 


1.5289e-2 


-2.096e-l 


3 


.005 


6.215 


-3.405e-4 


1.1219e-2 


-1.593e-l 


4 


.005 


6.455 


-2.343e-4 


8.0149e-3 


-1.182e-l 


5 


.005 


6.725 


-1.553e-4 


5.5331e-3 


-8.493e-2 


6 


.005 


7.025 


-9.92 le-5 


3.6904e-3 


-5.913e-2 


7 


.01 


7.37 


-5.083e-5 


1.9863e-3 


-3.344e-2 


8 


.035 


7.755 


-1.69 le-5 


7.0022e-4 


-1.252e-2 


9 


.155 


8.19 


-1.514e-6 


7.021 le-5 


-1.421e-3 


10 


.5 


8.675 


-6.080e-7 


3.1075e-5 


-6.9511e-4 


11 


.98 


9.22 


-3.832e-7 


2.0648e-5 


-4.867e-4 


12 


1.57 


9.84 


-1.168e-8 


6.6877e-7 


-1.669e-5 


13 


2.245 


10.53 


1.5193e-9 


-1.036e-7 


3.115e-6 


14 


3.005 


11.305 


-4.70e-10 


3.1455e-8 


-9.321e-7 


15 


3.855 


12.18 


-4.17e-10 


2.9655e-8 


-9.330e-7 


16 


4.81 


13.165 


-1.71e-10 


9.7399e-9 


-1.88e-7 


17 


5.885 


14.265 


3.798e-10 


-4.578e-8 


2.3450e-6 


18 


7.085 


15.505 


1.4311e-9 


-1.70e-7 


8.9045e-6 


19 


8.43 


16.895 


3.2174e-9 


-4.135e-7 


2.3613e-5 


20 


9.93 


18.46 


5.9566e-9 


-8.427e-7 


5.3146e-5 


21 


11.615 


20.215 


9.7702e-9 


-1.530e-6 


1.0702e-4 
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TABLE 15 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


2.4495e+0 


-1.104e+l 


3.1922e+l 


-5.942e+l 


6.9432e+l 


1 


2.0139e+0 


-9.343e+0 


2.7823e+l 


-5.329e+l 


6.4009e+l 


2 


1.6198e+0 


-7.76 le+0 


2.3864e+l 


-4.716e+l 


5.8369e+l 


3 


1.2752e+0 


-6.327e+0 


2.0135e+l 


-4.116e+l 


5.2596e+l 


4 


9.8193e-l 


-5.056C+0 


1.6691e+l 


-3.536e+l 


4.6763e+l 


5 


7.348 le-1 


-3.937e+0 


1.3521e+l 


-2.977e+l 


4.0872e+l 


6 


5.3395e-l 


-2.985e+0 


1.0692e+l 


-2.454e+l 


3.5082e+l 


7 


3.1748e-l 


-1.867e+0 


7.0416e+0 


-1.704e+l 


2.5757e+l 


8 


1.2640e-l 


-7.9357e-l 


3.2098e+0 


-8.392e+0 


1.3855e+l 


9 


1.6472e-2 


-1.2079e-l 


5.8335e-l 


-1.87 le+0 


3.9184e+0 


10 


8.9342e-3 


-7.2830e-2 


3.9193e-l 


-1.403e+0 


3.2817e+0 


11 


6.5912e-3 


-5.6606e-2 


3.2100e-l 


-1.212e+0 


2.9915e+0 


12 


2.3843e-4 


-2.1504e-3 


1.2706e-2 


-4.8814e-2 


1.1409e-l 


13 


-5.3993e-5 


5.9242e-4 


-4.2952e-3 


2.1380e-2 


-7.6762e-2 


14 


1.6489e-5 


-1.9865e-4 


1.7034e-3 


-9.777 le-3 


3.1192e-2 


15 


1.7543e-5 


-2.2598e-4 


2.1014e-3 


-1.3454e-2 


5.0992e-2 


16 


7.8184e-7 


2.0181e-5 


-3.0473e-4 


2.1525e-3 


-1.6839e-2 


17 


-6.7595e-5 


1.2191e-3 


-1.4467e-2 


1.1590e-l 


-6.3258e-l 


18 


-2.708 le-4 


5.2928e-3 


-6.9543e-2 


6.2408e-l 


-3.797e+0 


19 


-7.8867e-4 


1.7058e-2 


-2.4976e-l 


2.5099e+0 


-1.712e+l 


20 


-1.9670e-3 


4.7312e-2 


-7.7288e-l 


8.6887c+0 


-6.643e+l 


21 


-4.4011e-3 


1.1784e-l 


-2.147e+0 


2.6953e+l 


-2.304e+2 
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TABLE 15 (cont’d): EXPONENTIAL COEEHCIENTS, a = 0.948, b = 0.393 FOR 



SNR (dB) 


I 


J 


K 


0 


-4.861e+l 


1.9037e+l 


-4.423e+0 


1 


-4.595e+l 


1.8558e+l 


-4.660e+0 


2 


-4.308e+l 


1.8030e+l 


-4.960e+0 


3 


-4.001e+l 


1.7459e+l 


-5.342e+0 


4 


-3.676e+l 


1.6850e+l 


-5.827e+0 


5 


-3.334e+l 


1.6207e+l 


-6.444e+0 


6 


-2.983e+l 


1.5555e+l 


-7.227e+0 


7 


-2.337e+l 


1.3817e+l 


-8.087e+0 


8 


-1.418e+l 


■ 1.0848e+l 


-9.024e+0 


9 


-5.591e+0 


7.7044e+0 


-1.018e+l 


10 


-5.255e+0 


8.3106e+0 


-1.235e+l 


11 


-5.080e+0 


8.9005e+0 


-1.497e+l 


12 


-6.2182e-l 


5.5740e+0 


-1.664e+l 


13 


-2.8564e-l 


5.8552e+0 


-2.038e+l 


14 


-5.1125e-l 


6.8069e+0 


-2.537e+l 


15 


-5.5837e-l 


7.6168e+0 


-3.143e+l 


16 


-3.4763e-l 


8.0433e+0 


-3.858e+l 


17 


1.8330e+0 


4.3987e+0 


-4.363e+l 


18 


1.4609e+l 


-2.4846e+l 


-2.326e+l 


19 


7.5616e+l 


-1.873e+2 


1.5776e+2 


20 


3.3041e+2 


-9.572e+2 


1.1792e+3 


21 


1.2834e+3 


-4.199e+3 


6.069 le+3 
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TABLE 16: EXPONENTIAL COEFnCIENTS, a = 0.96043 b = 0.39783 FOR 



SNR (dB) 


Range 
Lower 
Limit XL 


Range 
Upper 
Limit xy 


A 


B 


C 


0 


.005 


5.715 


-7.7426e-4 


2.3443e-2 


-3.0587e-l 


1 


.005 


5.885 


-5.8026e-4 


1.8099e-2 


-2.4325e-l 


2 


.005 


6.075 


-4.2514e-4 


1.3692e-2 


-1.9001e-l 


3 


.005 


6.29 


-3.0242e-4 


1.0085e-2 


-1.4490e-l 


4 


.005 


6.535 


-2.0739e-4 


7.1837e-3 


-1.0721e-l 


5 


.005 


6.805 


-1.3815e-4 


4.9808e-3 


-7.7361e-2 


6 


.005 


7.115 


-8.7466e-5 


3.2954e-3 


-5.3479e-2 


7 


.01 


7.46 


-4.5137e-5 


1.7854e-3 


-3.0428e-2 


8 


.035 


7.85 


-1.5076e-5 


6.3191e-4 


-1.1431e-2 


9 


.155 


8.29 


-1.3689e-6 


6.4184e-5 


-1.3129e-3 


10 


.51 


8.785 


-5.4233e-7 


2.8055e-5 


-6.3508e-4 


11 


.995 


9.34 


-3.3557e-7 


1.8319e-5 


-4.3746e-4 


12 


1.595 


9.96 


-9.5417e-9 


5.5151e-7 


-1.3887e-5 


13 


2.28 


10.665 


1.3047e-9 


-9.0119e-8 


2.7462e-6 


14 


3.045 


11.450 


-4.066e-10 


2.7606e-8 


-8.2906e-7 


15 


3.91 


12.335 


-3.630e-10 


2.6193e-8 


-8.3554e-7 


16 


4.88 


13.330 


-1.531e-10 


8.9322e-9 


-1.8131e-7 


17 


5.965 


14.450 


3.247e-10 


-3.9778e-8 


2.0683e-6 


18 


7.18 


15.705 


1.2409e-9 


-1.4935e-7 


7.9300e-6 


19 


8.54 


17.115 


2.7948e-9 


-3.6393e-7 


2.1058e-5 


20 


1.0065 


18.695 


5.2091e-9 


-7.4668e-7 


4.7713e-5 


21 


11.77 


20.475 


8.5565e-9 


-1.3578e-6 


9.6223e-5 
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TABLE 16 (conl’d): EXPONENTIAL COEFFICIENTS, a = 0.96043 b = 0.39783 FOR 



SNR (dB) 


D 


E 


F 


G 


H 


0 


2.2505e+0 


-1.027e+l 


3.0060e+l 


-5.664e+l 


6.7007e+l 


1 


1.8439e+0 


-8.664e+0 


2.6135e+l 


-5.070e+l 


6.1684e+l 


2 


1.487 le+0 


-7.214e+0 


2.2459e+l 


-4.494e+l 


5.6311e+l 


3 


1.1740e+0 


-5.895e+0 


1.8988e+l 


-3.928e+l 


5.0801e+l 


4 


9.0203e-l 


-4.702e+0 


1.5715e+l 


-3.370e+l 


4.5126e+l 


5 


6.7728e-l 


-3.672e+0 


1.2761e+l 


-2.843e+l 


3.9496e+l 


6 


4.8909e-l 


-2.770e+0 


1.0046e+l 


-2.335e+l 


3.3813e+l 


7 


2.9237e-l 


-1.74 le+0 


6.6435e+0 


-1.627e+l 


2.4894e+l 


8 


1.1685e-l 


-7.4244e-l 


3.0390e+0 


-8.040e+0 


1.3431e+l 


9 


1.5385e-2 


-1.1400e-l 


5.5628e-l 


-1.803e+0 


3.8132e+0 


10 


8.2601e-3 


-6.8137e-2 


3.7103e-l 


-1.344e+0 


3.1805e+0 


11 


6.0012e-3 


-5.2212e-2 


2.9995e-l 


-1.147e+0 


2.8687640 


12 


1.9991e-4 


-1.8148e-3 


1.0772e-2 


-4.1387e-2 


9.54276-2 


13 


-4.8227e-5 


5.3596e-4 


-3.9345e-3 


1.9834e-2 


-7.2238e-2 


14 


1.4870e-5 


-1.8176e-4 


1.5822e-3 


-9.2194e-3 


2.984 le-2 


15 


1.5937e-5 


-2.0834e-4 


1.9667e-3 


-1.2785e-2 


4.9242e-2 


16 


1.0186e-6 


1.3693e-5 


-2.3094e-4 


1.6302e-3 


-1.40536-2 


17 


-6.0467e-5 


1.1053e-3 


-1.3287e-2 


1.0779e-l 


-5.95786-1 


18 


-2.4438e-4 


4.8394e-3 


-6.4419e-2 


5.8562e-l 


-3.609e+0 


19 


-7.1264e-4 


1.5617e-2 


-2.3165e-l 


2.3583e+0 


-1.630e+l 


20 


-1.7893e-3 


4.3607e-2 


-7.2179e-l 


8.2220640 


-6.369e+l 


21 


-4.0094e-3 


1.0877e-l 


-2.008e+0 


2.5545e+l 


-2.213e+2 
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TABLE 16 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.96043 b = 0.39783 FOR 



SNR (dB) 


I 


J 


K 


0 


-4.749e+l 


1.8824e+l 


-4.438e+0 


1 


-4.485e+l 


1.8342e+l 


-4.674e+0 


2 


-4.208e+l 


1.7826e+l 


-4.975e+0 


3 


-3.910e+l 


1.7267e+l 


-5.358e+0 


4 


-3.591e+l 


1.6660e+l 


-5.843e+0 


5 


-3.260e+l 


1.6031e+l 


-6.459e+0 


6 


-2.912e+l 


1.5373e+l 


-7.241e+0 


7 


-2.286e+l 


1.3669e+l 


-8,103e+0 


8 


-1.390e+l 


1.0743e+l 


-9.04 le+0 


9 


-5.492e+0 


7.6275e+0 


-1.020e+l 


10 


-5.150e+0 


8.2227e+0 


-1.237e+l 


11 


-4.936e+0 


8.7734e+0 


-1.498e+l 


12 


-5.7984e-l 


5.4749e+0 


-1.664e+l 


13 


-2.8179e-l 


5.7836e+0 


-2.039e+l 


14 


-4.9793e-l 


6.7184e+0 


-2.538e+l 


15 


-5.4500e-l 


7.5199e+0 


-3.144e+l 


16 


-3.4588e-l 


7.9523e+0 


”3.860c+l 


17 


1.7400e+0 


4.4387e+0 


-4.373e+l 


18 


1.4062e+l 


-2.412e+l 


-2.368e+l 


19 


7.2911e+l 


-1.829e+2 


1.5536e+2 


20 


3.2101e+2 


-9.423e+2 


1.1761e+3 


21 


1.2489e+3 


-4.140e+3 


6.0649e+3 
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V. COMPARISON OF THE APPROXIMATIONS 
aMax[\I\,\Q,\} + bMin{\l\,\Q\} AND al/l + b\Q\ OF A 
TARGET TEST CELL 



Signals detected by the envelope detection approximations and given by (1) and 
(2) have significantly different probability density functions. This fact is true for both cases 
of noise and noise plus target. A graphical comparison of x^ and x^ with the density 



function of the true envelope detector x 



-J. 



I^ + Q^ is shown in Figure 30 for the case 




Figure 30: PDF for x, x^ and Xg, No Target, /i = 1 

in which no target is present. These curves are generated by numerical convolution using 
the program in Appendix D. The density of Xg is a result of the convolution of (46) and (48). 
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The density of Xg is obtained by convolving equations (83) and (84). Finally, the density 
function of x is obtained from [Ref. 8] 



Pj]T ^2 = xe ^u(x) . (99) 

The coded sections for each of these are identified in the program. The multiplying 
coefficients of both approximations are arbitrarily set at a = 1.0 and b = 0.25. 

The results of the numerical convolutions shown in Figure 30 are verified against 
experimental results. Data sets generated for / and Q are sorted into histograms, shown in 
Figure 31. The star symbols “ * ” represent the true envelope detector jc, while the thinner 




Value of Trials 

Figure 31: Histograms of x, x^ and Xg, No Target, n = 1 

bars represent Xg and the wider bars represent Xg. The program in Appendix E is used to 
generate these histograms. 
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As shown in Figure 30, the envelope detector approximation models the true 
envelope detector the best. The approximation Xg is significantly shifted from the true 
detector. These curves are used to compute the function Py ^ in (7) and (8), for the Pd and 
the PpA respectively. As explained in Section n.C, Py „ is an n-fold convolution of these 
curves, for a given n reference cells. 

The disparity between Xg and x for the function Py „ is enhanced as « is increased. 
This can be seen in Figures 32 ad 33, where each of x, Xg and Xg has been convolved 




Figure 32: PDF for x, Xg and Xg, No Target, n = 4 

for n = 4 and n = 16. These n-fold convolutions are generated using the program at 
Appendix F. The approximation Xg generally has the same mean as the true detector x. The 

mean of the approximation Xg, however, is increasingly offset from that of the true detector 
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as n increases. For n= 1, the mean of is offset from x by about 0.5 dB in SNR. For n = 
4, this offset is about 1.25 dB, whereas for n = 16, the offset is approximately 4 dB. 




Figure 33: PDF for x, x^ and x^, No Target, n = 16 



Such differences ultimately create significant differences in the Pj) and the Pp^ GO 
CFAR performance, as compared to the values obtained for a true envelope detector. As an 
example, the experimental Pj) curves for each of jc, x^ and x^ are plotted in Figure 34. These 

are generated using the Monte Carlo simulation program at Appendix G for a = 1.0 and b 
= 0.25 and n = 4. This results of this simulation program are validated against similar 
simulations performed in previous research [Ref. 3, Ref. 4]. The thresholds for each 

detector are chosen so as to yield a PpA - lO"^. The threshold for is selected from 
available plots [Ref. 1] as T. = 5.537. The thresholds for x and x are determined from 
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experimental trials using the program at Appendix H. They are set at = 3.5 and 
r„ = 4.2. The Pd curves show that Xg provides the most degraded detection sensitivity 




Figure 34: Probability of Detection Curves for x, x^ and Xg, a = 1.0, b = 0.25, n = 4. 

compared to the true envelope detector. The approximation Xg provides a significantly 
better Pd performance in comparison. As n is increased, this advantage becomes even more 
pronoimced. 
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VI. CONCLUSION 



A. SUMMARY 

The main contribution of this thesis is to determine if the probability of detection for 
a GO CFAR radar processor using the envelope detection approximation 
= aMax { |/1, \Q \ } + bMin { |/|, \Q \ } can be expressed analytically in a closed form. 
The solution to this problem depends solely on whether the probability density function of 
a radar range ceU can be derived in a closed form for both of the cases where a target is 
present and where no target is present. It is shown in Section in that this cannot be done 
without the use of either an error function approximation or an exponential curve-fitting of 
the density function. 

The closed form solution based on an error function approximation is very complex. 
As stated in Section HI. A.6.a, it is comprised of more than 1 100 first order exponential 
terms. For the closed form evaluation of the each of these terms needs to be integrated 

in closed form. This approach is not practical. Therefore, the best option is to evaluate the 
Pp numerically once the error fvmction approximation has been applied to Although this 

is stUl very complex and tedious to program, it has the advantage of being one of the most 
accurate methods of computing the Pp analytically, as long as the residual error in the error 
function approximation is kept to a minimum. Another advantage is that this approach is 
valid for any value of SNR and multiplying coefficients a and b. It can therefore also be 
used to compute the Ppy^. 

The exponential curve-fitting approach involves the curve-fitting of a number of 
curves generated numerically for a given range of signal-to-noise ratios. The exponential 
coefficients resulting from these curve-fits are then summarized in reference tables. This 
has the advantage of enabling a rapid implementation of p^ within the Pp expression. Then 
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the Pj) can be evaluated numerically, as per the error function approximation approach. The 
main disadvantage of this method is its lack of versatility. It is dependent upon the 
multiplying coefficients a and b of the envelope approximation. Furthermore it can only be 
solved for whole number values of SNR (i.e. 0, 1,2, ), upto the SNR required for the 

Pd=1- 

An additional envelope detection approximation Xg = a\I\+b\Q\ is considered in 
Section rV. It is offered as an alternative to x^. Since it is a simplified version of Xg, its 
density function Pg (x) is expressed relatively simply in a closed form solution. It is also 
expressed in the form of exponential curve-fits (as done for Xg), with the required 
coefficients summarized in reference tables. 

Finally, the two envelope detection approximations x^ and Xg are compared to each 
other and to the true envelope detector x. This is done in terms of the density functions for 
non-target range cells and in terms of the Pj) performance. In both instances, it is shown 

that the approximation x^ delivers a performance which is significantly better than that 
provided by Xg. Ultimately, the choice of envelope detection approximation to be used is 
to be based on Pj) performance and the ease of implementation. Although Xg performs 

better in terms of detection capabilities, Xg is simpler to implement since the Max and Min 
operations do not have to be performed. 

B. FUTURE RESEARCH RECOMMENDATIONS 

The possibilities for future research on GO CFAR performance using the envelope 
detection approximation Xg = aMax { |/| , |(2I } + bMin {\I\,\Q\] are qmte attractive from 
a mathematical point of view. The ultimate goal of such research would be to provide the 
radar design and analysis community with a simple method to compute the Pj) and the Pp^ 
theoretically. This solution would have to be independent of the signal-to-noise ratio in the 
test cell (or reference cell) and independent of the multiplying coefficients a and b. 
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Simplified analytical solutions for the Pj) and the Pp^ are invariably contingent upon 

the ability to integrate expressions of the form of equations (77) and (78) in closed form. 
This presents a challenging mathematical problem. Solving it would certainly prove to be 
very useful to anyone using this particular approximation Xg for an envelope detector, 
whether it is for a GO CFAR processor or for some other use. 

Although less glamorous in its possibilities, the ability to approximate the error 
function in a single segment — «> to while maintaining a very small residual error would 
be very useful. This would greatly reduce the complexity of the closed form probability 
density function of a test cell as presented in Section in.6. This could be done by 
appropriately modulating either an inverse tangent function or a hyperbolic tangent 
function. A closed form solution would then be required for the integral of this function 
multiplied by an exponential of second order. 

Finally, the generation of the Pj) curves using the analytical solutions presented in this 
thesis would be very useful. This could be done for both types of envelope detection 
approximations Xg and Xg, for all the cases listed in Table 1. The number of reference cells 

used should be n = 1, 2, 4, 8, 16 and 32. These represent the majority of the situations 
encountered in GO CFAR processing. 
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APPENDIX A - MATLAB PROGRAM FOR 
POLYNOMIAL CURVE-FITTING OF ERROR 

FUNCTION 



%polynomial curve-fitting of error function 

x=.01:.01;.5; 

y=erf(x); 

pl=polyfit(x,y,2); 

f=polyval(pl,x); 

e=abs(f-y); 

plot(x,y,x,f,’;’),grid 

xlabel(‘x’) 

ylabel(‘erf(x)’) 

%subplot(212) 

%plot(x,e),grid 
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APPENDIX B - C PROGRAM FOR NUMERICAL 
CONVOLUTION FOR OVER A RANGE OF SNR 

/^probability density functions of various envelope detectors 
for a GO CFAR */ 

#include <stdio.h> 

#include <math.h> 

#define pi (4.0*atan(1.0)) 

main(){ 

FILE *Pcell, *PcellArea; 

int idelay,c,cl,c2,x,dataLength,s; 

float n; 

double pdfCellTgt[16000],pdfSample,pMax[ 16000], 
pMin[ 1 6000], phi, sum, 

u,snr,ss,a,b,m 1 ,m2,e 1 ,e2,e3,e4,ef 1 ,ef2,ef3,ef4; 

dataLength= 1 6000; 
pdfSample=l .0/200.0; 
a=1.0; 
b=0.25; 

if ((Pcell=fopen(“Pyn2.dat”,”w”))=NULL) { 
printfC'????”); 
exit(l); } 

for (s=0;s<=29;s++){ 

snr=pow( 10.0,((float)s/10.0)); 
ss=sqrt(2*snr); 

for (c=0;c<dataLength;c++) pdfCellTgt[c]=0.0; 

for (phi=0;phi<=9*pi/40;phi+=pi/40){ 
ml=ss*cos(phi); 
m2=ss*sin(phi); 

for (c=0;c<dataLength;c++){ 
u=(double)c*pdfS ample; 

/*envelope detection approximation aMax{III,IQI}+bMin{III,IQI} */ 
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el=exp(-pow((u/a-ml),2.0)/2.0)+exp(-pow((u/a+ml),2.0)/2.0); 
e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
e3=exp(-pow((u/b-m l),2.0)/2.0)+exp(-pow((u/b+m 1 ),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
efl=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2=erf((u/a-m 1 )/sqrt(2.0))+erf((u/a+m 1 )/sqrt(2.0)); 
eB=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2.0-erf((u/b-ml)/sqrt(2.0))-erf((u/b+ml)/sqrt(2.0)); 
pMax[c]=(el*efl + e2*ef2)/(2.0*a*sqrt(2.0*pi)); 
pMin[c]=(e3*ef3 + e4*ef4)/(2.0*b*sqrt(2.0*pi)); 

/♦envelope detection approximation alll + bIQI */ 

/* el=exp(-pow((u/a-ml),2.0)/2.0)+exp(-pow((u/a+ml),2.0)/2.0); 

e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin [c] =e4/(b *sqrt(2 .0*pi)) ; 
pMax[c]=el/(a*sqrt(2.0*pi)); */ 

/♦envelope detector sqrt(III^2 + IQI'^2) ♦/ 

/♦ pdfCellTgt[c]=u*exp(-pow(u,2.0)/2); ♦/ 

) 

idelay=0; 

for (c=0;c<dataLength;++c) { 
idelay+=l; 
sum=0.0; 

for (cl=l;cl<=idelay;++cl) 
sum+=(pMax[cl ] ♦pMin[idelay-c 1 + 1 ] ); 
pdfCellTgt[c]+=(sum^pdfS ample); 

} 



j ^****************gp^ Qp ph^se IJQQP********************* f 

for (x=0;x<dataLength;x++){ 
pdfCellTgt[x]=pdfCellTgt[x]/10.0; 
fprintf(Pcell,”%e\n”,pdfCellTgt[x]); 

} 



J y*****gJ>J£) QP ^QQp**************************^ 

fclose(Pcell); 

^4:************5f:*************3H************** :^:4:*3)c**** 2ND MAIN */ 
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APPENDIX C - MATLAB PROGRAM FOR 
EXPONENTIAL CURVE-FITTING OF FOR A 

RANGE OF SNR 



%curve-fitting of target test cell pdf curves over 
%a range of SNR values (in dB) 

load el 

c=el; 

clg 

clear p 
hold 

order=10; 

SNR=21; 

sampleSize= 1/200; 
for q=l:SNR+l, 
rl(q)=l; 

while c(rl(q),q)<.00001, 
rl(q)=rl(q)+l; 
end 

r2(q)=rl(q)+100; 
while c(r2(q),q)>.00001, 
r2(q)=r2(q)+l; 
end 

x=[r 1 (q)- 1 :r2(q)- 1] ’ .*sampleSize; 

p(q,:)=polyfit(x,log(c(rl(q):r2(q),q)),order); 

f=exp(polyval(p(q,:),x)); 

plot(x,c(rl(q):r2(q),q),x,f,’--’,x,abs(f-c(rl(q):r2(q),q))) 

end 

hold 
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APPENDIX D - C PROGRAM FOR WITH NO 

TARGET 



/*probability density functions of various envelope detectors 
for a GO CFAR for no target case */ 

. #include <stdio.h> 

#include <math.h> 

#define pi (4.0*atan(1.0)) 

mainQl 

FILE *Pcell, *PcellArea; 

int idelay,c,cl,c2,x,dataLength,s; 

float n; 

double pdfCellTgt[16000],pdfSample,pMax[16000], 
pMin[ 16000],phi,sum, 

u,snr,ss,a,b,m 1 ,m2,e 1 ,e2,e3,e4,ef 1 ,ef2,ef3,ef4; 

d ataLength= 1 6000; 
pdfS ample= 1 .0/200.0; 
a=1.0; 
b=0.25; 

if ((Pcell=fopen(“Pyn2.dat”,”w”))==NULL){ 
printf(“????”); 
exit(l); } 



ss=0.0; 



for (c=0;c<dataLength;c-H-) pdfCellTgt[c]=0.0; 



ml=ss; 

m2=ss; 



for (c=0;c<dataLength;c++){ 
u=(double)c*pdfS ample; 

/♦envelope detection approximation aMax{III,IQI}+bMin{III,IQI} */ 

e l=exp(-pow((u/a-m l),2.0)/2.0)+exp(-pow((u/a+m l),2.0)/2.0); 
e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
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e3=exp(-pow((u/b-ml),2.0)/2.0)+exp(-pow((u/b+ml),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
efl=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2=erf((u/a-m 1 )/sqrt(2.0))+erf((u/a+m 1 )/sqrt(2.0)); 
ef3=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2.0-erf((uA>-ni 1 )/sqrt(2.0))-erf((u/b+m l)/sqrt(2.0)); 
pMax[c]=(el*efl + e2*ef2)/(2.0*a*sqrt(2.0*pi)); 
pMin[c]=(e3*ef3 + e4*ef4)/(2.0*b*sqrt(2.0*pi)); 

/♦envelope detection approximation alll + bIQI */ 

/* el=exp(-pow((u/a-ml),2.0)/2.0)+exp(-pow((u/a+ml),2.0)/2.0); 

e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin[c]=e4/(b*sqrt(2.0*pi)); 
pMax[c]=el/(a*sqrt(2.0*pi)); */ 

/♦envelope detector sqrt(III^2 + IQI^'2) ♦/ 

/♦ pdfCellTgt[c]=u^exp(-pow(u,2.0)/2); ♦/ 

} 

idelay=0; 

for (c=0;c<dataLength;++c){ 
idelay+=l; 
sum=0.0; 

for (c 1 = 1 ;c 1 <=idelay;++c 1 ) 
sum+=(pMax[c 1 ] ♦pMin[idelay-cl + 1 ] ); 
pdfCeUTgt[c]+=(sum^pdfS ample); 

} 

for (x=0;x<dataLength;x++){ 
pdfCellTgt[x]=pdfCellTgt[x]/10.0; 
fprintf(Pcell,”%e\n”,pdfCeUTgt[x]); 

} 

fclose(Pcell); 

J END MAIN */ 
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APPENDIX E - MATLAB PROGRAM COMPARISON 
OF ENVELOPE DETECTION APPROXIMATION 



METHODS BY HISTOGRAM 

%comparison of various envelope detectors 
%for GO CFAR by histogram 

clear 

clg 

b=.25; 

a=l; 

N=30000; 

q=randn(l,N); 

i=randn(l,N); 

Q=abs(q); 

I=abs(i); 

z=sqrt(i.''2 + q.''2); 
z 1 =a*max(I,Q)+b*min(I,Q); 
z2=a*I + b*Q; 

[m,n]=hist(z,200); 

[m 1 ,n 1 ] =his t(z 1 ,200) ; 

[m2,n2]=hist(z2, 1 00); 

bar(n,m./((l/48)*N),’*’) 

hold 

bar(n 1 ,m 1 7(( 1/48)*N)) 
bar(n2,m2./((l/24)*N)) 
hold 
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APPENDIX F - C PROGRAM FOR n-FOLD 



CONVOLUTION FOR Py „ 

/*n-fold convolution of reference cell PDF for a GO-CFAR */ 

#include <stdio.h> 

#include <math.h> 

#define pi (4.0*atan(1.0)) 

main(){ 

FILE *Pcell, *PcellArea; 
int idelay,c,cl,c2,x,dataLength; 
float s,n; 

double pdfCellTgt[ 1 6000] ,pdfS ample.phaseS ample,pMax[ 1 6000], 
p2[16000],pl[16000],p3[16000],p4[16000],p5[16000], 
pMin[ 1 6000],phi,theta,sumPhase[ 16000],area,sum, 
u,snr,ss,a,b,m 1 ,m2,e 1 ,e2,e3,e4,ef 1 ,ef2,eB,ef4; 

d ataLength= 1 6000; 
pdfSample= 1 .0/200.0; 
a=1.0; 
b=.25; 

if ((PceU=fopenC‘pyn.dat”,”w”))=NULL) { 
printf(“????”); 
exit(l); } 



ss=0.0; 

ml=ss; 

m2=ss; 



for (c=0;c<dataLength;c++){ 
u=(double)c*pdfSample; 

/♦envelope detection approximation aMax{III,IQI}+bMin{III,IQI} */ 

/* el =exp(-pow( (u/a-m 1 ),2 .0)/2 .0)+exp(-pow((u/a+m 1 ) ,2 .0)/2.0) ; 

e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
e3=exp(-pow((u/b-ml),2.0)/2.0)+exp(-pow((u/b+ml),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
efl=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2=erf((u/a-m 1 )/sqrt(2.0))+erf((u/a+m l)/sqrt(2.0)); 
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ef3=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2 .0-erf ((u/b-m 1 )/sqrt(2 ,0))-erf ((u/b+m 1 )/sqrt(2 .0)); 
pMax[c]=(el*efl + e2*ef2)/(2.0*a*sqrt(2.0*pi)); 
pMin[c]=(e3*eO + e4*ef4)/(2.0*b*sqrt(2.0*pi)); */ 

/*enveIope detection approximation alll + bIQI */ 

/* el=exp(-pow((u/a*ml),2.0)/2.0)+exp(-pow((u/a+ml),2.0)/2.0); 

e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin[c]=e4/(b*sqrt(2.0*pi)); 
pMax[c]=el/(a*sqrt(2.0*pi)); */ 

/*true envelope detector sqrt(III'^2 + IQI'^2) */ 

pi [c]=u*exp(-pow(u,2.0)/2); 



} 

/*do not perform this convolution for true envelope detector */ 
/* idelay=0; 

for (c=0;c<dataLength;++c){ 
pl[c]=0.0; 
idelay+=l; 
sum=0.0; 

for (cl=l;cl<=idelay;++cl) 
sum+=(pMax [c 1 ] * pMin [idelay-c 1+ 1 ] ) ; 
p 1 [c]=sum*pdfSample; 

} *! 

/*convolve up to here for n=l */ 

idelay=0; 

for (c=0;c<dataLength;++c){ 
p2[c]=0.0; 
idelay+=l; 
sum=0.0; 

for (cl=l;cl<=idelay;-H-cl) 
sum+=(p 1 [c 1 ] *p 1 [idelay-c 1 + 1 ] ); 
p2[c]=sum*pdfSample; 

} 

/*convolve up to here for n=2 */ 
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idelay=0; 

for (c=0;c<dataLength;++c){ 
p3[c]=0.0; 
idelay+=l; 
sum=0.0; 

for (cl=l;cl<=idelay;++cl) 
sum+=(p2[cl]*p2[idelay-cl+l]); 
p3[c]=sum*pdfSample; 

} 

/*convolve up to here for n=4 */ 
idelay=0; 

for (c=0;c<dataLength;++c){ 
p4[c]=0.0; 
idelay+=l; 
sum=0.0; 

for (cl=l;cl<=idelay;++cl) 
sum+=(p3[c 1 ]*p3 [idelay-c 1 + 1 ]); 
p4[c]=sum*pdfSample; 

/* fprintf(PceU,”%e\n”.p4[c]); */ 

} 

/♦convolve up to here for n=8 */ 
idelay=0; 

for (c=0;c<dataLength;++c){ 
p5[c]=0.0; 
idelay+=l; 
sum=0.0; 

for (c 1 = 1 ;c 1 <=idelay ;++c 1 ) 
sum+=(p4[cl]*p4[idelay-cl+l]); 
p5[c]=sum*pdfSample; 
fprintf(PceU,”%e\n”,p5[c]); 

} 

/♦convolve up to here for n=16 ♦/ 
fclose(Pcell); 

J END MAIN */ 
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APPENDIX G - MATLAB PROGRAM FOR MONTE 



CARLO SIMULATIONS OF Pp 

%GO CFAR Monte Carlo simulation for probability of detection 
clear 

T=3.5; 

Tl=4.2; 

T2=5.537; 

a=l; 

b=.25; 

n=4; 

S=10000; 

PHI=pi/4; 

phi=0:PHI/20:PHI; 

snr=0:25; 

w=2*n; 

for B=l:length(phi), 
for C=l:length(snr), 

SNR=10'^(snr(C)/10); 

A=sqrt(2*SNR); 

d=0; 

dl=0; 

d2=0; 

X=zeros(l,S); 

I=A*cos(phi(B)) + randn(l,S); 

Q=A*sin(phi(B)) + randn(l,S); 

X=sqrt(I.''2 + Q.^2); 

Xl=a*max(abs(I),abs(Q)) + b*min(abs(I),abs(Q)); 

X2=a*abs(I) + b*abs(Q); 

i=randn(S,w); 
q=randn(S,w); 
x=sqrt(i.^2 + i.^2); 

xl=a*max(abs(i),abs(q)) + b*min(abs(i),abs(q)); 
x2=a*abs(i) + b*abs(q); 

ld=sum(x(:,l:n)’); 
lg=sum(x(: ,n+ 1 : w) ’ ); 
z=max(ld,lg); 
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Vt=T*z/n; 



ldl=sum(xl(:,l:n)’); 

lgl=sum(xl(;,n+l;w)’); 

zl=max(ldl,lgl)J 

Vtl=Tl*zl/n; 

Id2=sum(x2(:,l:n)’); 
Ig2=sum(x2(:,n+1 :w)’); 
z2=max(ld2,lg2); 
Vt2=T2*z2/n; 

d=length(find((X-Vt)>0)); 

Pd(B,C)=d/S; 

dl=length(find((Xl-Vtl)>0)); 

Pdl(B,C)=dl/S; 

d2=Iength(find((X2-Vt2)>0)); 

Pd2(B,C)=d2/S; 

end 

end 

PD=mean(Pd); 

PDl=mean(Pdl); 

PD2=mean(Pd2); 

save PD PD PDl PD2 
quit 



APPENDIX H - MATLAB PROGRAM FOR MONTE 



CARLO SIMULATIONS OF PpA 

%GO CFAR Monte Carlo simulation for probability of 

%false alarm 

clear 

T=3.5; 

Tl=4.2; 

a=l; 

b=25; 

n=4; 

S=100000; 

w=2*n; 

d=0; 

dl=0; 

X=zeros(l,S); 

I=randn(l,S); 

Q=randn(l,S); 

X=sqrt(I.'^2 + Q.'^2); 

Xl=a*max(abs(I),abs(Q)) + b*min(abs(I),abs(Q)); 

i=randn(S,w); 
q=randn(S,w); 
x=sqrt(i.'^2 + i.^2); 

xl=a*max(abs(i),abs(q)) + b*min(abs(i),abs(q)); 

ld=sum(x(:,l;n)’); 

lg=sum(x(:,n+l:w)’); 

z=max(ld,lg); 

Vt=T*z/n; 

ldl=sum(xl(:,l:n)’); 

lgl=sum(xl(:,n+l;w)’); 

zl=max(ldl,lgl); 

Vtl=Tl*zl/n; 

d=length(fmd((X-Vt)>0)); 

Pfa=d/S; 

dl=length(fmd((Xl-Vtl)>0)); 
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Pfal=dl/S; 



save PFA Pfa Pfal 
quit 
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